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We report on measurements of dielectron (e + e _ ) production in Au+Au collisions at a center- 
of-mass energy of 200 GeV per nucleon-nucleon pair using the STAR detector at RHIC. Sys¬ 
tematic measurements of the dielectron yield as a function of transverse momentum (px) and 
collision centrality show an enhancement compared to a cocktail simulation of hadronic sources 
in the low invariant-mass region (M ee < lGeV/c 2 ). This enhancement cannot be reproduced 
by the p-meson vacuum spectral function. In minimum-bias collisions, in the invariant-mass 
range of 0.30 — 0.76GeV/c 2 , integrated over the full pt acceptance, the enhancement factor is 
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1.76 ± 0.06 (stat.) ± 0.26 (sys.) ± 0.29 (cocktail). The enhancement factor exhibits weak centrality 
and pt dependence in STAR’S accessible kinematic regions, while the excess yield in this invariant- 
mass region as a function of the number of participating nucleons follows a power-law shape with 
a power of 1.44 ± 0.10. Models that assume an in-medium broadening of the p meson spectral 
function consistently describe the observed excess in these measurements. Additionally, we report 
on measurements of oj and (j >-meson production through their e + e _ decay channel. These mea¬ 
surements show good agreement with Tsallis Blast-Wave model predictions as well as, in the case 
of the 0-meson, results through its K + K~ decay channel. In the intermediate invariant-mass re¬ 
gion (1.1< M ee < 3 GeV/c 2 ), we investigate the spectral shapes from different collision centralities. 
Physics implications for possible in-medium modification of charmed hadron production and other 
physics sources are discussed. 

PACS numbers: 25.75.Cj, 25.75.Dw 


I. INTRODUCTION 

A major scientific goal of the ultra-relativistic 
heavy-ion program is to study Quantum Chromo- 
Dynamics (QCD) matter at high temperature and 
density in the laboratory. Previous measurements 
at Relativistic Heavy Ion Collider (RHIC) have es¬ 
tablished the formation of a strongly-coupled Quark 
Gluon Plasma (sQGP) in high-energy heavy-ion col¬ 
lisions [lj. Throughout the evolution of the hot, 
dense, and strongly interacting system, electromag¬ 
netic probes are produced and escape with little 
interaction. Thus, they provide direct information 
about the various stages of the system’s evolution. 

Following convention, the dilepton invariant-mass 
spectrum is typically divided into three ranges: the 
low mass region - LMR (Mu < M^), the interme¬ 
diate mass region - IMR (M^ < Mu < Mj/,p) and 
the high mass region - HMR (Mu > Mj/^). As 
will be described next, distinctively different phys¬ 
ical processes contribute or even dominate within 
these particular ranges. 

The initial hard perturbative QCD process, Drell- 
Yan production, ( qq —> l + l~) can produce high-mass 
dileptons and is expected to be an important mecha¬ 
nism in the HMR [2|. Moreover, initial hard scatter¬ 
ing processes can allow for bremsstrahlung emission 
of virtual photons which convert into low invariant 
mass, high transverse momentum (py) dielectrons 
(“internal conversion”). These dileptons, in prin¬ 
ciple, are calculable within the perturbative QCD 
framework. 

The colliding participant system is expected to 
quickly reach the partonic sQGP phase where dilep¬ 
tons can be produced through electromagnetic ra¬ 
diation via parton-parton scatterings. Theoretical 
calculations indicate that at top RHIC energy, QGP 
thermal dilepton production will become a domi¬ 
nant source in the IMR while thermal dileptons with 
higher masses originate from earlier stages Q ■ This 
suggests that investigating the thermal dilepton pro¬ 
duction as a function of Mu & Pt allows for prob¬ 
ing the medium properties at different stages of the 


system’s space-time evolution. Measuring thermal 
dilepton collective flow and polarization can reveal 
information about the relevant degrees of freedom 
which may relate to deconfinement and the equi¬ 
librium of the strongly interacting matter created 
in heavy-ion collisions Thermal radiation can 

produce real photons as well as virtual photons that 
decay to dileptons. Comparative analysis of distri¬ 
butions for these dileptons with respect to those pro¬ 
duced in initial hard scattering, can shed light on di¬ 
rect real photon production from the QGP medium. 

When the system expands, cools down, and en¬ 
ters into the hadronic phase, dileptons are produced 
via multiple hadron-hadron scattering by coupling 
to vector mesons (p, w, <p, etc.). They are ex¬ 
pected to dominate the LMR and their mass spectra 
may be related to the chiral symmetry restoration 
in the medium Si- Theoretical calculations sug¬ 
gest that the vector meson spectral functions will 
undergo modifications in a hot and dense hadronic 
medium, which may be connected to the restoration 
of chiral symmetry. Two scenarios have been pro¬ 
posed for the change of vector meson spectral func¬ 
tions when chiral symmetry is restored: a shift of 
the pole mass 0 and/or a broadening of the mass 
spectrum jTl]. Measurements of the diclectron con¬ 
tinuum in the low mass region will help expose the 
vector meson production mechanisms, and hence the 
chiral properties of the medium in heavy-ion colli¬ 
sions. 

Finally, when all particles decouple from the sys¬ 
tem, long-lived hadrons (7r°, p, DD } etc.) can decay 
into lepton pairs and are measured by the detector 
system. Their contributions can be calculated based 
on the measured or predicted invariant yields of the 
respective parent particles, and incorporated in the 
so-called hadron cocktail. 

Dilepton measurements in heavy-ion collisions 
have been pursued for decades from relatively low 
ener gies to relativistic and ultrarelativistic ener¬ 
gies The CERES measurements of e + e“ 

spectra in Pb+Au collisions at the Super Proton 
Synchrotron (SPS) showed an enhancement in the 
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mass region below ^700 MeV/c 2 with respect to the 
hadron cocktail that included the vacuum line shape 
for the p meson [il|- High-statistics measurements 
from the NA60 experiment at ^/snn = 17.2 GeV sug¬ 
gested that this enhancement is consistent with in¬ 
medium broadening of the p-meson spectral func¬ 
tion rather than a drop of its pole mass E m- 
l2lj . Strikingly, after removal of correlated charm 
contributions, the NA60 collaboration also observed 
that the slope parameters of the dimuon transverse 
mass (tot) spectrum showed a roughly linear in¬ 
crease with dimuon invariant mass below the (f>- 
meson mass, followed by a sudden decline at higher 
masses. This observation provided a first indication 
of thermal leptons from a partonic source fl 6 ] . 

Thermal radiation of dileptons is expected to 
be significantly enhanced due to a well-developed 
partonic phase in the heavy-ion collision systems 
created at RHIC. The PHENIX collaboration has 
measured dielectrons at mid-rapidity in Au+Au 
collisions within its detector acceptance. For 
minimum-bias collisions, in the mass region be¬ 
tween 150 and 700 MeV/c 2 an enhancement of 
4.7 ± 0.4 (stat.) ± 1.5 (syst.) ±0.9 (model) has been 
reported by the PHENIX collaboration [17]. Several 
theoretical calculations, which have successfully ex¬ 
plained the SPS data flU (3 and the STAR data [22} 
previously, were unable to reproduce the magnitude 
of the low mass dielectron enhancement observed 
by PHENIX through expected vector meson con¬ 
tributions in the hadronic medium. The PHENIX 
measured IMR yields are consistent with the charm 
contribution from p + p scaled with the number of 
binary collisions. However, within the limits of the 
data precision and our present understanding of the 
modification of charmed hadron production in Au 
+ Au collisions, no conclusive evidence for thermal 
radiation can be inferred from this measurement. 
A detailed dilepton program to investigate the in¬ 
medium chiral and thermal properties is one of the 
main focuses of future heavy-ion projects from Schw- 
erionen Synchrotron (SIS) energies up to LHC ener¬ 
gies. 

In this paper, we report on detailed measurements 
of dielectron production in Au + Au collisions at 
^/snn = 200 GeV with the Solenoidal Tracker At 
RHIC (STAR) experiment. The data used in this 
analysis were recorded during the RHIC runs in 2010 
and 2011. The barrel Time-Of-Flight (TOF) detec¬ 
tor system was completed before these runs thus sig¬ 
nificantly improving the electron identification over 
a wide momentum range in STAR’S Time Projection 
Chamber (TPC). 

The paper is organized as follows: Section II de¬ 
scribes the experimental setup and the data sets 
used in this analysis. Section III explains in detail 


the analysis techniques including electron identifica¬ 
tion, dielectron invariant mass reconstruction, back¬ 
ground subtraction, detector acceptance efficiency 
correction, and systematic uncertainties. Section IV 
presents our results on dielectron production yields 
within the STAR detector acceptance and a compar¬ 
ison to the hadron cocktails. Results are compared 
with theoretical calculations of in-medium modified 
vector meson line shapes as well as QGP thermal 
radiation contributions. Systematic studies on the 
centrality and px dependence of the dielectron yields 
are presented. Our results and conclusions are sum¬ 
marized in Section V. 


II. EXPERIMENTAL SETUP 

The data used in this analysis were collected by 
the STAR detector [23| . The major detector subsys¬ 
tems used in this analysis are the TPC, the TOF, 
and two trigger subsystems: the Vertex Position De¬ 
tectors (VPDs) and the Zero Degree Calorimeters 
(ZDCs). 

A. Time Projection Chamber 

The TPC [24] is the main tracking detector and 
consists of a 4.2 m long solenoidal cylinder concen¬ 
tric with the beam pipe. It is operated in a uniform 
0.5 Tesla magnetic field parallel to the beam direc¬ 
tion (defined as z direction in STAR). The inner and 
outer radii of the active volume are 0.5 and 2.0 m, 
respectively. It covers the full azimuth and a pseu¬ 
dorapidity range of |ry| < 2 for the inner radius and 
1 77 1 < 1 for the outer radius. The TPC has 45 read¬ 
out layers allowing measurements of charged particle 
momenta with a resolution of ~ 1% at px~ 1 GeV/c 
for tracks originating from the collision vertices. It 
is also used for particle identification (PID) via the 
ionization energy loss ( dE/dx) in the TPC gas with 
a mean dE/dx resolution of about 7% . 

B. Time Of Flight System 

The TOF system consists of the Barrel TOF 
(BTOF) detector covering the TPC cylinder and the 
VPDs located in the forward pseudorapidity regions. 
The latter provide the common start time. BTOF 
detector utilizes the multi-gap resistive plate cham¬ 
ber technology (25|. It covers the full azimuth and 
| r ]| < 0.9. The VPD detector has two parts, sitting 
along the beam pipe on both sides of the STAR de¬ 
tector at ±5.7 m from the center. The detectors 
cover a pseudorapidity range of 4.4< \rj\ <5.1 [26| . 
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The time stamps recorded by the VPD and the 
BTOF detectors are used to calculate the particle 
time-of-flight (tof). The tof is further combined with 
the track length and momentum, both measured by 
the TPC, to identify charged particles. The timing 
resolution of the TOF system, including the start 
timing resolution in Au+Au 200 GeV collisions, is 
less than 100 ps. 


C. Trigger Definitions 

The minimum bias trigger in Au+Au 200 GeV col¬ 
lisions for the 2010 and 2011 runs was defined as 
a coincidence between the two VPDs and an on¬ 
line collision-vertex cut in order to select collision 
events that took place near the center of the detec¬ 
tor. The central trigger in the 2010 Au+Au col¬ 
lisions includes the ZDC detectors, located on both 
sides of the STAR detector at approximately ±18 m. 
This trigger requires a small signal in the ZDC de¬ 
tectors in combination with a large hit multiplicity 
in the BTOF and corresponds to the top 10% of the 
total hadronic cross section. 


III. ANALYSIS TECHNIQUE 
A. Event Selection and Centrality Definition 

Events used in this analysis were required to have 
a reconstructed collision vertex (primary vertex) 
within 30 cm of the TPC center along the beam di¬ 
rection in order to ensure uniform TPC acceptance. 
To suppress the chance of selecting the wrong ver¬ 
tex from different bunch-crossing collisions and to 
ensure that the selected event indeed fired the trig¬ 
ger, the difference between event vertex z-coordinate 
K TPC and the K VPD calculated from the VPD tim¬ 
ing was required to be within 3 cm. These se¬ 
lection criteria yield 240M (year 2010) and 490M 
(year 2011) 0-80% minimum-bias triggered events 
and 220M (year 2010) central triggered (0-10%) 
Au+Au events at ^/snn = 200 GeV. The results re¬ 
ported in this paper are from the combined year 2010 
and year 2011 data. 

Centrality in Au+Au 200 GeV collisions was de¬ 
fined using the uncorrected charged particle multi¬ 
plicity dN/dr] within \r)\ < 0.5. The dN/dr) distri¬ 
bution was then compared to a Monte Carlo (MC) 
Glauber calculation in order to delineate the cen¬ 
trality bins. Furthermore, the dependence of dN/dr] 
on the collision vertex position V z and the beam lu¬ 
minosity has been included in order to take accep¬ 
tance and efficiency changes on the measured dN/dr) 
into account. The measured uncorrected dN/dr) 



FIG. 1. Upper panel: Uncorrected multiplicity 
distribution measured within \r)\ < 0.5 and \V Z \ < 5 cm. 
The solid curve depicts the multiplicity distribution from 
a MC Glauber simulation. Bottom panel: Ratio between 
MC and data. 

distribution from Au+Au 200 GeV minimum-bias 
events collected in year 2010 is shown in Fig. [I] 
The dN/dr) distributions are from the V z region of 
—5 < V z < 5 cm and extrapolated to a zero ZDC- 
coincidence rate, so as to correct for the detector 
acceptance and efficiency dependence on the V z and 
luminosity. The measured distribution matches the 
MC Glauber calculation well for dN/dr) > 100. In 
the lower multiplicity region, the VPD trigger be¬ 
comes less efficient. The bottom panel shows the ra¬ 
tio between MC and measured data. The centrality 
bins are defined according to the MC Glauber dis¬ 
tribution in order to determine the centrality cut on 
the measured dN/dr). To obtain the real minimum- 
bias sample, events in the low multiplicity region 
have been weighted with the ratio shown in Fig. |T] 
(bottom panel) to account for the VPD inefficiency. 

The average number of participants (-Vp art ), and 
number of binary collisions {N\,- ln ) from MC Glauber 
simulations of Au+Au at ^/snn = 200 GeV are 
listed in Table Q] 


B. Track Selection 

Electron candidate tracks used in this analysis 
were required to satisfy the following selection crite¬ 
ria: 

• the number of fit points in the TPC (nHits- 
Fits) should be greater than 20 (out of a max- 
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TABLE I. Summary of centrality bins, average number 
of participants (Apart), and number of binary collisions 
(Abin) from MC Glauber simulation of Au+Au at ^/snn 
= 200 GeV. The errors indicate uncertainties from the 
MC Glauber calculations. 


Centrality 

(Apart) 

(Abin) 

0-10% 

325.5 + 3.7 

941.2 ± 26.3 

10-40% 

174.1 ± 10.0 

391.4 + 30.3 

40-80% 

41.8 + 7.9 

56.6 + 13.7 

0-80% 

126.7 + 7.7 

291.9 + 20.5 


imum of 45) to ensure good momentum reso¬ 
lution; 

• the ratio of the number of fit points over the 
number of possible points should be greater 
than 0.52 in order to avoid track splitting in 
the TPC; 

• the distance of closest approach (DCA) to the 
primary vertex should be less than 1 cm in 
order to reduce contributions from secondary 
decays; 

• the number of points used for calculating 
( dE/dx) (nHitsdEdx) should be greater than 
15 to ensure good dE/dx resolution; 

• the track should match to a valid TOF hit with 
the projected position within TOF’s sensitive 
readout volume. 


Au + Au \ s NN = 200 GeV 



(b) 


2x10' 1 1 2 

momentum (GeV/c) 


FIG. 2. (Color online) Upper panel: normalized dE/dx 
(no e ) vs. momentum ( p ) distributions for all charged 
particles. Bottom panel: na e vs. p distributions after 
applying the TOF velocity cut |l//3 — 1| < 0.025. 


C. Electron identification 


Electrons (including positrons if not specified) 
were identified based on a combination of the TPC 
and TOF detectors. The electron identification pro¬ 
cedure has been described in [27[ . In low multiplicity 
collisions, electrons can be cleanly separated from 
hadrons by requiring a TOF velocity cut and us¬ 
ing the TPC truncated mean ionization energy loss 
dE/dx dependence on particle momentum. How¬ 
ever, the situation becomes more complicated in 
high multiplicity Au+Au collisions. The normalized 
dE/dx is defined as follows: 


na e 


In ((dE/dx) m /{dE/dx) f) 


R, 


dE / dx 


(i) 


where (} m and (} th represent measured and theoret¬ 
ical values, respectively, and RdE/dx is the experi¬ 
mental dE/dx resolution. The na e vs. p distribu¬ 
tions for the 2010 data are shown in Fig. [2] The 
upper panel shows the distribution for all charged 
particles, the bottom panel shows the distribution 


after applying the TOF velocity selection 11//3 — 11 < 
0.025, which accepts about 95% of the electrons 
based on the TOF timing resolution. Despite the 
TOF velocity selection, there are still some slow 
hadrons that contribute to the electron band in this 
distribution. The source of these remaining slow 
hadrons is described in the following paragraphs. 

The na e distribution for TPC tracks with matched 
TOF hits are shown in Fig. [2] For most cases where 
TOF hits are correctly associated with the charged 
particle tracks, one would have a meaningful parti¬ 
cle velocity measurement that can then be used for 
particle identification. There are also many TOF 
hits from electrons that originate from photon con¬ 
versions in the material between the TPC and TOF 
sensitive detector volumes. Since photons do not 
leave a trace in the TPC, these TOF hits can be 
randomly associated with TPC tracks especially in 
high-multiplicity events. 

The inverted particle velocity (1//3) measured by 
the TOF (time) and the TPC (path length) versus 
the particle momentum ( p ) measured by the TPC 
is shown in the upper panel of Fig. [3] for all TPC- 
TOF associations in Au+Au collisions at Asnn = 
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200 GeV. The band below 1//3 = 1 depicts the asso¬ 
ciations between conversion electron TOF hits and 
random TPC charged tracks. The bottom panel of 
Fig. ED shows the 1//3 distributions in the momen¬ 
tum range 0.2 < p < 0.25 GeV/c for three centrality 
classes. The three distributions are normalized to 
the pion peak. One can see that with increasing mul¬ 
tiplicity that the fake association fraction increases 
substantially. These random associations were fur¬ 
ther confirmed using Monte Carlo (MC) Geant [28 ] 
simulations. 

As mentioned before, the TOF-based velocity of 
particles depends on the time-of-flight measurements 
from the TOF detector and the track length deter¬ 
mined by the TPC. For particles from secondary 
vertex decays ( e.g ., 7r, AT, and p from Kg, A, and 12 
decays), the track length and time-of-flight measure¬ 
ments have some offset, which leads to uncertainties 
when calculating the velocity. 

Consequently, the applied particle velocity cut 
cannot remove the random association of charged 
hadron tracks with TOF signals and the particles 
from secondary vertex decays. Such hadrons are 
mostly at momentum of 400 MeV/c or above where 
the hadron dE/dx bands cross the electron band. 
These hadrons remain in the dE/dx vs. p distribu¬ 
tion in the lower panel of Fig. [2] and introduce an 
additional hadron background in the sample of se¬ 
lected electron candidates in the region where the 
electron dE/dx band crosses with hadrons (mostly 
kaons and protons). The dashed black lines in the 
lower panel of Fig. [2] depict the dE/dx cuts to se¬ 
lect the single electron candidates in this analysis. 
Finally, distributions of the number of selected elec¬ 
tron candidates are shown in upper panel of Fig. [I] 
their raw pt spectra are also shown in the lower 
panel of Fig. 0] 


D. Electron Purity and hadron contamination 

The ncr e vs. p distribution after the TOF velocity 
selection has been shown in Fig. [2] We have per¬ 
formed a multi-component fit to the na e distribu¬ 
tion for individual momentum slices in order to de¬ 
compose the yields of each particle species, and thus 
derive the electron purity and hadron contamination 
for a certain na e cut. The na e distribution for elec¬ 
trons is assumed to be Gaussian, with its position 
and shape determined by selecting conversion elec¬ 
trons using an invariant-mass reconstruction. The 
positions and shapes of the na e distributions for pi- 
ons, kaons, and protons were determined by selecting 
pure samples of these particles with particle masses 
calculated from the TOF detector. Figure [5] shows 
the respective m 2 thresholds for the pure hadron 



Momentum (GeV/c) 



FIG. 3. (Color online) Upper panel: 1//3 vs. particle 
momentum. Solid line depicts a prediction for those as¬ 
sociations where TOF hits were triggered by conversion 
electrons while matched randomly with TPC charged 
tracks. Bottom panel: 1//3 projection in the momen¬ 
tum bin 0.2 < p < 0.25 GeV/c for three centrality bins, 
normalized to the pion peak region. 

samples. The positions and shapes of all compo¬ 
nents are kept fixed during the fits, leaving only the 
individual yields as free parameters to fit the na e dis¬ 
tribution slices in Fig. [5] An example of the fit result 
for the momentum bin of 0.68 < p < 0.72 GeV/c is 
shown in Fig. [6] The black dotted curve at high na e 
region depicts a small contribution of tracks from 
merged pions in the TPC. The (dE/dx) value of 
these tracks are twice that of normal pion tracks, 
thus its position and shape is predictable from the 
pion na e distribution. For completeness, we have in¬ 
cluded this contribution in the fit although it is well 
separated from the electron peak. 

The multi-component fits describe the full distri- 












nc e 

FIG. 6. Upper panel: na e distribution of clean ir,K,p 
samples selected using TOF masses. Dashed lines are 
Gaussian fits to these distributions allowing the ex¬ 
trapolation of the tail regions where contamination be¬ 
comes apparent. Bottom panel: An example of multi- 
component fit to the ncr e distribution for the momen¬ 
tum bin 0.68 < p < 0.72 GeV/c in 200 GeV Au+Au 
minimum bias collisions. The two vertical lines indicate 
the selected electron range. 


FIG. 4. Upper panel: The distributions of number of 
electron candidates. Bottom panel: Raw pr spectra of 
the electron candidates. 



m 2 (GeV/c 2 ) 


butions well in the regions where the slow hadron 
peaks can be separated from the electron peaks. In 
the region where kaons and protons start to overlap 
with electrons, we use the hadron yields from neigh¬ 
boring momentum bins with clean particle identifi¬ 
cation to interpolate the expected hadron yield. The 
systematic uncertainties on the electron purity in 
these overlapping bins were estimated by compar¬ 
ing the yields to the results from the free paramter 
fit, which take the hadron yield as free parame¬ 
ter. Figure [7] shows the electron purity for the 
candidate samples used in the minimum-bias and 
central collisions. As expected, hadron contamina¬ 
tion increases from peripheral to central collisions. 
The electron purity integrated over the region of 
0.2<p T <2.0 GeV/c is (94.6±1.9)% and (92.1±2.0)% 
for 0-80% minimum bias and 0-10% central Au+Au 
collisions, respectively. The impact of hadron con¬ 
tamination on the dielectron spectra will be further 
discussed in Section IIII HI 


FIG. 5. Charged particle m 2 distribution from TOF 
measurements in 200 GeV Au+Au minimum bias colli¬ 
sions. The shaded areas are the respective m 2 thresholds 
used for selecting high purity 7r, K,p samples. 


E. Electron pairing and background 
subtraction 

For each individual event, all electron and 
positron candidates within the STAR acceptance of 
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momentum (GeV/c) 


FIG. 7. (Color online) Estimated electron purity vs. mo¬ 
mentum in 200 GeV Au + Au collisions. Gray areas in¬ 
dicate the momentum regions where na e of kaons and 
protons cross with that of electrons resulting in large 
uncertainties in those ranges. 

Pt> 0.2 GeV/c and \r]\ <1 are combined to generate 

the inclusive unlike-sign pair (iV_|_) invariant-mass 

distribution. Despite slight acceptance differences 
between the TPC and the TOF, the collision vertex 
distribution along beam direction ( z ) will provide fi¬ 
nite acceptance and efficiency for charged tracks up 
to |?y| <1. Therefore we used \iq\ <1 in this analy¬ 
sis and the dependence of efficiency and acceptance 
along t) has been corrected in the final spectra. In 
Fig. [ 8 j a two-dimensional distribution in invariant- 

mass and transverse-momentum (M ee , pr) of N _)_ 

pairs is shown in the STAR acceptance with \y ee \ <1 
(electron pair rapidity) and the background sub¬ 
tracted. Vector meson signals (oj, 4>, and J/^) are 
fairly easy to recognize after the background sub¬ 
traction. All distributions shown in this paper are 
calculated within the same STAR acceptance includ¬ 
ing \y ee \ <1 unless specified otherwise. 

In this analysis the signal ( S ) is defined as the 
e + e _ pairs that originate from pair production 
sources such as 7 r°, 77 , 77 ', p, w, <j>, J/t/>, 7 * decays, as 
well as correlated charm hadron decay. Background 
sources that contribute to the inclusive unlike-sign 
pair distributions include: 

• Combinatorial background pairs from two un¬ 
correlated electrons. 

• Correlated background pairs. For instance, in 
the case of Dalitz decays followed by a conver¬ 
sion of the decay photon(e.< 7 ., 7r° —> e + e~ r y, 
then 7 Z —> e + e~Z*), the electron from the 
Dalitz decay and the positron from the conver¬ 
sion are not completely uncorrelated as they 
originate from the same source. Another sig¬ 



M ee (GeV/c 2 ) 


FIG. 8. Two-dimensional (M ee , Pt) distribution of 
unlike-sign e + e _ pairs with background subtraction from 
200 GeV Au + Au minimum bias collisions in the STAR 
acceptance (pT>0.2GeV/c, |?7|<1, and \y ee \<l). 

nificant contribution is the electron pairs from 
same-jet fragmentation or back-to-back di-jet 
fragmentation. This source may become more 
significant at high mass or px ■ 

Contributions from uncorrelated and correlated 
background pairs are thoroughly studied and eval¬ 
uated using like-sign pairs, N++ and N —, con¬ 
structed from the same event. It has been demon¬ 
strated that when the e + and e~ are produced in sta¬ 
tistically independent pairs, the geometric mean of 

the like-sign pairs 2^/N++ x N _fully describes the 

background in the inclusive unlike-sign pair distribu¬ 
tion N .|_ 0 - In this analysis, we consistently used 

the like-sign distribution 2^/N ++ x N _to esti¬ 

mate or normalize the background distribution. The 
mixed-event unlike-sign distribution B _| was con¬ 
structed to estimate the combinatorial background 
and was used for a better statistical background 
estimation wherever the correlated background is 
negligible or the mixed-event unlike-sign distribu¬ 
tion agrees with the same-event like-sign distribution 
2^/N ++ x N _ Mixed-event like-sign pair distri¬ 

butions B ++ , were also constructed to verify 
the applicable kinematic region for the mixed-event 
technique as well as to define the normalization fac¬ 
tor for the mixed-event unlike-sign distribution. 

A sizeable component of the correlated electron 
pairs, that is not considered as part of the final sig¬ 
nal distribution, originates from photon conversions 
in the detector material. Details of the conversion 
electron removal will be discussed in subsection III- 
E.l. 

Hadron contamination in the selected elec¬ 
tron/positron sample due to particle misidentifica- 
tion may result in some residual contributions to the 
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final signal distribution. Most of these are from res¬ 
onance decays. The high purity of the electron sam¬ 
ple in this analysis allows us to demonstrate that the 
residual contribution due to hadron contamination 
in the final distribution is negligible. Such details 
will be discussed in the Section III-H. 


1. Photon conversion removal 

Background pairs from photon conversion were re¬ 
moved from the sample using the </>v angle selection 
method. This method is similar to that used by the 
PHENIX collaboration [l7] and relies on the kine¬ 
matics of the pair production process. The opening 
angle between the two conversion electrons should 
be zero, and the electron tracks are bent only in the 
plane perpendicular to the magnetic held direction, 
which for the STAR experiment is parallel to the 
beam direction (z). Unit-vector definitions used for 
the construction of the (f>y angle were taken from [|17] 
as: 


u = 


P+ +P- 
I P+ +P -1 


,v=p+ X 


U) = UXV,U) C = UXZ 


COS </>v = u> ■ w c 


( 2 ) 


where p± are momentum vectors of tracks, and 
z is the magnetic held direction. 

For pairs that originate from photon conversions 
(fry should be zero. It has no preferred orientation for 
combinatorial pairs, and only very weak dependence 
for e + e _ pairs from hadron decays. The electron 
pair mass versus (j>y for conversion electron pairs 
from the full Geant simulation of the STAR detec¬ 
tor [l^] is shown in the upper panel of Fig. [9] The 
populated bands at different mass positions depict 
the conversion electron pairs from different detec¬ 
tor materials. The reconstructed masses are shifted 
from zero as the electrons are assumed to originate 
from the primary vertex during the hnal track recon¬ 
struction. As a result, the three main bands from 
low to high masses correspond to the conversions 
from the beam pipe (at a radius r ~ 4 cm), inner 
cone support structure (r ~ 20 cm), and TPC inner 
held cage (IFC) (r ~ 46 cm). In order to remove 
these conversion pairs, we dehne a mass-dependent 
( t>v selection which is shown as the red line in the 
upper panel of Fig. [9l We estimated that more than 
95% conversion pairs are removed by this selection 
criterion. 

The signal pair invariant mass spectra before and 
after this photon conversion cut are shown in the 
bottom panel of Fig. [9] their difference is shown 



Au + Au \fs NN = 200 GeV (MC) 

— mass vs b cut 


0.6 0.8 

4> v (rad) 



FIG. 9. Upper panel: (f>y vs. mass distributions for pho¬ 
ton conversion electron pairs form the full Geant MC 
simulation. The solid red line depicts the mass depen¬ 
dent (f>y cut that was used to remove these conversion 
pairs. Bottom panel: photon conversion contribution in 
200 GeV Au+Au minimum bias collisions. The insert 
plot shows the structures from the beam pipe, the sup¬ 
porting bars of the inner cone, and the TPC inner field 
cage. 


as the filled histogram. Like-sign background sub¬ 
traction was used to obtain these distributions. Al¬ 
most all conversions appear in the mass region below 
0.1GeV/c 2 . 

The cut removing the photon conversion pairs was 
applied only in the very low mass region (M ee < 
0.2 GeV/c 2 ). The effect of the cut on the mixed- 
event distribution normalization is negligible as that 
determination is done at a much higher mass region. 


2. Event mixing 

The event mixing technique was used to reproduce 
the combinatorial background with improved statis¬ 
tical precision. In order to make the mixed-event 
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distributions close to that from real events, we have 
only selected events with similar properties for the 
mixed-event calculation. The full sample is divided 
into different pools according to the following event 
level properties: multiplicity, vertex position, event 
plane angle, and magnetic field direction. The sort¬ 
ing by event multiplicity and vertex position ensures 
electrons are mixed between events with similar de¬ 
tector acceptance and efficiency. This technique 
has been widely used in many other STAR analy¬ 
ses for the reconstruction of the combinatorial back¬ 
grounds [29j . The small signal-to-background ratio 
requires a very good understanding of the mixed- 
event distribution in the dielectron analysis. Its de¬ 
pendence on the event pool division for event-plane 
angle and magnetic field direction were studied in 
detail and are presented here. 

Elliptic flow measurements (30] in 200 GeV 
Au+Au collisions have shown that the momentum 
phase space distribution of particles produced in the 
event is approximately elliptical. Therefore, we only 
mix events with similar event-plane direction to en¬ 
sure the events have similar momentum phase space 
alignment, and further guaranteed by the multiplic¬ 
ity assortment to ensure the events have similar mo¬ 
mentum phase space distributions. The event plane 
was reconstructed with a conventional method using 
tracks in the TPC (0.1 < pr < 2 GeV/c and \tj\ < 1) 
in order to obtain the second-order event-plane an¬ 
gle + [30], In Fig. [TT3 results of a study using min¬ 
imum bias Au+Au collisions in which mixed-event 
unlike-sign and like-sign distributions are compared 
using different numbers of event pools in event-plane 
angle are shown. The figure illustrates several sce¬ 
narios from 1 up to 24 event pools. The dashed 
lines at ±0.5% corresponds to a 100% change in the 
yield where the signal-to-background ratio is 1/200. 
This study shows the importance of doing the divi¬ 
sion in event-plane angle in order to avoid distortion 
of the mixed-event distributions. The distortion is 
quite clear in the low mass region (< 1 GeV/c 2 ) 
and not negligible in the intermediate mass region 
(1—3 GeV/c 2 ). The differences become negligible 
when the number of event pools is 12 in 200 GeV 
Au±Au minimum-bias collisions, comparable to the 
TPC 2 nd order event plane resolution. 

A similar study of the centrality dependence for 
background distributions was carried out. As a re¬ 
sult, to ensure the minimal difference in all centrality 
bins studied, we choose 24 event pools in the event- 
plane angle in our analysis. 

The data samples used in this analysis were taken 
under two different magnetic field configurations of 
similar magnitude but opposite direction. The ac¬ 
ceptance for oppositely charged tracks in the two 
magnetic field configurations is not exactly the same 



FIG. 10. (Color online) Difference in mixed-event dis¬ 
tributions when using different number of event pools in 
the reconstructed event-plane angle. The left plot shows 
differences of unlike-sign distributions with different di¬ 
visions and the right plot shows differences of ratios of 
unlike-sign to like-sign distributions. 


due to a slight offset of the beam line with respect to 
the center of STAR detector system. Only electrons 
from events with the same magnetic field configu¬ 
ration were mixed when constructing total mixed- 
event distributions. 

The final number of event pools used in track 
multiplicity, vertex position, event plane angle, and 
magnetic field configuration isl6xl0x24x2 for 
this analysis of the 200 GeV minimum-bias Au±Au 
data. 

The statistics in the mixed-event distributions de¬ 
pend on the number of events chosen for the calcu¬ 
lation. In practice, however, the calculation can be 
done to sufficient precision for every event pool with 
a sizable number of event pools under the limitations 
of the number of events during the calculation. The 
differences between mixed-event distributions with 
different number of events in the buffers are shown 
in Fig. HU We observe no distortions beyond statis¬ 
tics in our calculation using a buffer of 50 events per 
event pool. With this choice, the statistical uncer¬ 
tainties in the mixed-event background are negligible 
compared to the same event distributions. 


3. Mixed-event normalization 

The unlike-sign and like-sign pair distributions in 

the same event (N^ _, N ++ / _), and in the mixed- 

event ( B _|_, B ++ / _) were constructed in two di¬ 

mensions ( M ee , pt)- The mixed-event unlike-sign 
distribution (£?_]_) provides the shape of the un¬ 

correlated combinatorial background, with an over¬ 
all normalization factor determined separately. The 
normalization factor was evaluated from the like- 
sign pair distribution using the technique described 
in Ref. E3- This technique is susceptible to a sys- 
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FIG. 11. (Color online) Ratios between mixed-event dis¬ 
tributions using different number of events to be mixed 
with in event buffers. 


tematic bias if correlated pairs exist in the the like- 
sign sample. Therefore, the kinematic region used 
to evaluate the normalization is carefully selected 
where such correlated backgrounds are negligible. 

The procedure to obtain the normalized combina¬ 
torial background !Tj°™ b is described in Ref. and 
also shown in the following Eq. [3] 


= 


A- = 


Jn.r. N ++ (M,p T )dMdp T 
Jn.r. B ++ (M,p T )dMdp T 
Jn.r. N — {M,p T )dMdp T 


B 


++ 


Jn.r. 
/»00 

JO 


px)dMdpx 

A + B ++ ( M , p T )dM dpT 


B 1 ' 


B 


comb / 


/»oo 

/ A_B _ (M,pT)dMdpT 

Jo 

2 y/B. 


norm . ^norm 


+ _ (M, pt) j™ B+ _ dMdpT 


B + ~(M,pt) 


( 3 ) 


N.R. denotes the integral calculated in a certain 
kinematic region, i. e. the normalization region. Ta¬ 
ble un lists the total like-sign pairs in the normaliza¬ 
tion region for each centrality class and the corre¬ 
sponding statistical uncertainties of the normaliza¬ 
tion factors. 


TABLE II. Total Like-sign pairs in the normalization re¬ 
gion (N.R.) in each centrality class and the correspond¬ 
ing statistical uncertainties of the normalization factors. 


Centrality 

Like-sign pairs in N.R. 

Statistical un. 

0-80% 

4.2 x 10° 

4.8 x 10 -4 

0-10% 

8.9 x 10° 

3.3 x 10" 4 

10-40% 

2.3 x 10° 

6.5 x 10" 4 

40-80% 

8.0 x 10 5 

1.1 x 10 ;i 


The residual difference between same-event like- 
sign iV++._and the normalized mixed-event 



M ee (GeV/c 2 ) 4 



FIG. 12. (Color online) Upper panel: Residual differ¬ 
ences between same-event and mixed-events like-sign dis¬ 
tributions divided by its standard deviation. The black 
box indicates the default normalization region. Bottom 
panel: Residual difference distributions for all the entries 
in different mass regions. 


as a function of M ee and px is shown in the 
upper panel of Fig. [T3] The difference is normal¬ 
ized by the expected statistical error in each kine¬ 
matic bin. The residual difference distributions for 
all entries in different mass regions are shown in the 
bottom panel of Fig. [T2] In the black box in the up¬ 
per panel of Fig. [12l the normalized residuals follow 
the statistical fluctuation. We then chose this area 
1 < M ee < 2 GeV/c 2 as the normalization region in 
our analysis. The systematic uncertainty introduced 
by the selected normalization region was studied by 
varying the selection as will be discussed in more 
detail in Section III-H. 

In Fig. 1131 the raw mass distributions of mixed- 
event like-sign and unlike-sign pairs in the full px 
region are plotted together with the same event dis¬ 
tributions. To further investigate any residual differ¬ 
ences between these distributions, the ratios between 
them are plotted in Fig. [TT] Panels (a-c) show that 
in the normalization region the residuals are negligi¬ 
ble. The slight increasing trend in the higher mass 
region can be attributed to the possible jet-related 
correlated background E3- This will be discussed 
further in Section III-E.5. The px and centrality 
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FIG. 13. (Color online) Raw pair mass distributions for 
200 GeV Au + Au collisions. The mixed-event unlike- 
sign and like-sign distributions are normalized in the 
mass region from 1.0 to 2.0GeV/c 2 . 


dependence of the inclusive unlike-sign and the nor¬ 
malized mixed-event mass distributions are shown 
in Fig. [15] 


4. Like-sign and unlike-sign acceptance difference 
correction 

The like-sign distribution is taken as the best esti¬ 
mate for the background in the inclusive unlike-sign 
distribution. However, the acceptances for like-sign 
and unlike-sign pairs differ in the STAR detector 
due to the magnetic field. The observed candidate 
e + and e~ tracks <f> versus px are shown in Fig. 1161 
The empty strips along the <f> direction are due to 
the TPC read-out sector boundaries. These accep¬ 
tance boundaries and local inefficiencies or accep¬ 
tance holes in the active detecting area will results 
in different acceptances for like-sign and unlike-sign 
pairs. We used the mixed-event technique to calcu¬ 
late these acceptance differences. 

The correction factor for the acceptance difference 
between like-sign and unlike-sign pairs is obtained as 
a ratio of the like-sign and unlike-sign distribution 
from mixed-event. The ratio was calculated in each 
(M ee , px) bin, and the corresponding correction ap¬ 
plied in this 2D plane. The geometric mean from the 

two like-sign charge combinations H—|-,-describes 

the background in the unlike-sign -)— combinations 
in total pairs in spite of any detecting efficiency E3- 
When calculating the combined like-sign pair in each 
kinematic bin, we use both the geometric mean and 
the direct sum of +-1- and-pairs in the calcula¬ 

tion to estimate the impact of potentially different 
detecting efficiencies for positive and negative tracks, 
shown in Eqs. []] and [5] 


N±± (M,pr) = 2 \f JV++(M,pr) • N —(M,px) 


_ B+-{M 1 p T ) _ 

2 • yjB ++ (M, pt) ■ B — (M,pr) 


(4) 


and 


N^(M,p T ) = a[N ++ (M,p T ) + IY__(M,p T )] 


_ B+-(M,pr) _ 

b • [ B ++ (M,p T ) + B — (M,p T )\ 


_ /o°° 2 • y/N ++ (M, PT ) ■ N__(M,p T )dMd PT 
f^°[N ++ (M,p T ) + N —( M,p T )]dMdp T 
b _ f 0 °° 2 • \JB ++ (M, p T ) • B — (M,'p^)dMdpT 

/ 0 °°[H ++ (M,p T ) + B — (M,p T )]dMdp T 


(5) 


where N ++ , 1V__, B ++ , and B _denote the dis¬ 
tributions of like-sign (++) and (-) from the 

same event and mixed-event calculation, respec¬ 
tively. !?_| denotes the unlike-sign distribution 

from mixed-event calculations. N±± denotes the 
acceptance-corrected like-sign background distribu¬ 


tion. 

In Fig. 1171 the ratio of mixed-event unlike-sign 
and like-sign distributions is shown as a function of 
the pair mass integrated over px. The structures ob¬ 
served in the ratio at low mass are caused by local 
inefficiencies and acceptance holes. This ratio has a 
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FIG. 14. Panel (a), (b), (c): Ratios between same event and mixed-event like-sign distributions. 
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FIG. 15. (Color online) pr (left panel) and centrality (right panel) dependence of the same-event inclusive unlike-sign 
distributions (histograms) and the normalized unlike-sign mixed-event background distributions (red lines). 
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FIG. 16. (Color online) (j> vs. pr for all negative (left 
panel) and positive (right panel) tracks from a single 
magnetic field configuration. The blank areas are due 
to the TPC sector boundaries, which shows the different 
acceptances between positive and negative tracks partic¬ 
ularly in the low pr due to the magnetic field. 


acceptance hole due to the merging is significantly 
smaller than the estimate that is used. As a result, 
we conclude that effects due to track merging in the 
TPC are negligible. 

Signal loss can also occur when two TPC tracks 
point to the same TOF read-out cell (size 6x3 cm 2 at 
a typical radius of about 215 cm). The TOF match¬ 
ing algorithm removes any TPC-TOF association in 
this situation since it cannot resolve the timing of 
two close hits. To evaluate such losses, pairs are arti¬ 
ficially removed for which the TPC tracks pointed to 
neighboring TOF cells, thereby increasing the accep¬ 
tance hole by a factor of about 9. The impact on the 
final acceptance correction factor is ~0.05% and lim¬ 
ited to two particular mass regions (^0.35 GeV/c 2 
due to unlike-sign pairs, ~0.1 GeV/c 2 due to like- 
sign pairs). 



FIG. 17. Acceptance correction factor for unlike-sign and 
like-sign pair difference from 200 GeV Au+Au minimum- 
bias collisions. 

dependence on the pair px and a correction is ap¬ 
plied to the like-sign distributions in the 2D (M ee , 
pr) plane. 

There are additional inefficiencies from merging 
effects that are different for like-sign and unlike- 
sign pairs in a magnetic field. These inefficiencies 
can originate from TPC-track merging or TOF-hit 
merging. We use two-particle correlations to study 
this acceptance loss due to the TPC-track merging. 
We calculate the Ary and A (f> correlations of like- 
sign and unlike-sign pairs in both same and mixed 
events. As a conservative estimation, we artificially 
remove a significant amount of the detection area 
near (Ar] 1 A<p)= (0,0), and correct the background- 
subtracted spectra with the cut efficiency which was 
estimated by the mixed events. The difference in the 
final mass spectrum was <1%. The actual TPC hit 
resolution is around 1 mm, for which the expected 


5. Correlated background 

In this analysis, the like-sign distribution is used 
as the best estimate of the background in the in¬ 
clusive unlike-sign distribution. The properly nor¬ 
malized mixed-event unlike-sign distributions were 
taken as the combinatorial background contribution. 
The difference between the like-sign and the mixed- 
event unlike-sign was used to understand the corre¬ 
lated background contributions. 

The ratio of the acceptance corrected like-sign to 
the mixed-event unlike-sign distributions is shown 
in Fig. [TU In the low-mass region (<1 GeV/c 2 ), 
the difference is due to the cross-pair contributions 
such as 7 r° —> e + e~ r y, followed by 7 Z —> e + e~Z*. 
In the intermediate and high-mass regions, the like- 
sign and mixed-event distributions generally agree 
within our current precision, but also show a trend 
of an increasing excess with increasing mass. This 
trend is expected to be mostly due to back-to-back 
jet correlations. 

We use a data-driven method to estimate the cor¬ 
related background contribution. We fit the ratio in 
Fig. HU in the mass region above 1 GeV/c 2 with two 
different empirical functions: a second order poly¬ 
nomial and an exponential function. The small dif¬ 
ference from unity in these fits is assigned as resid¬ 
ual correlated background. We use the 68.3% con¬ 
fidence limits from the fit Eq. [ 6 ] (indicated by the 
dashed lines in the figure) as the systematic uncer¬ 
tainty on the correlated background. The lower limit 
of this uncertainty is consistent with unity, indicat¬ 
ing that the like-sign background is consistent with 
the mixed-event unlike-sign background. 

R{M) = 1 + e (M " a)/6 


( 6 ) 
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FIG. 18. Ratio of the same-event like-sign to the mixed- 
event unlike-sign distributions. The gray area indicates 
the normalization region. The solid and dashed lines de¬ 
pict an empirical fit to the distribution in the mass region 
of 1—4GeV/c 2 and the fit uncertainties, respectively. 


This residual background has been studied in dif¬ 
ferent pt and different centrality bins. The accep¬ 
tance correction factors, which are estimated via the 
ratio between unlike-sign and like-sign mixed-event 
distributions are shown in Fig. [121 In Fig- Ell the ra¬ 
tios of acceptance corrected likc-sign backgrounds to 


mixed-event unlike-sign distributions are also shown 
for various px and centrality selections. The accep¬ 
tance correction factor shows a slight centrality de¬ 
pendence as the number of electron candidates is 
different in each centrality. On the other hand, it 
shows a strong px dependence due to the varying 
track curvatures in the magnetic field for tracks as 
a function of px- At sufficiently high px, tracks are 
nearly straight, and the acceptance of like-sign and 
unlike-sign pairs is expected to be similar. A data- 
driven procedure was used to estimate the correlated 
background in each px and centrality bin. 


6. Signal extraction 

In this analysis, the dielectron signal for invariant 
masses of M ee < 1.0 GeV/c 2 is obtained by subtract¬ 
ing the same-event like-sign background from the in¬ 
clusive unlike-sign distribution. In the higher mass 
region, we first subtract the combinatorial back¬ 
ground using the mixed-event unlike-sign pairs for 
better statistical precision. The residual correlated 
background is evaluated by the data-driven method 
described in the previous subsection and subtracted 
together with the combinatorial background. The 
signal extraction evaluated over the entire invariant 
mass region reported here is described as follows: 


N + _(M,pt) — N±±(M,pt) forM<M th 

N + -(M, P t) - B™ffi h (M,p T ) x [1 + r(M,px)] for M > M th 


where r(M,px) is the correlated background contri¬ 
bution normalized to the mixed-event combinatorial 
background and M t h is 1.0 GeV/c 2 in our default 
calculations. We vary this transition mass point be¬ 
tween 1.0—2.0 GeV/c 2 and find the difference in the 
final mass spectrum to be negligible (<0.05%). 


The raw signal invariant mass spectrum, 

S'-!_(M,px), for 200 GeV Au+Au minimum-bias 

collisions obtained by applying Eq. Q] is shown in 
the top panel of Fig. [20] along with the inclusive 
unlike-sign and background distributions. The 
bottom panel shows the signal-to-background ratio 
(S/B) in p+p [13 and Au+Au collisions. For the 
latter, the S/B at M ee = 0.5 GeV/c 2 is about 
1/200 in minimum-bias and 1/250 in 0-10% central 
collisions. 


F. Efficiency and Acceptance Correction 

The raw dielectron signal yields must be corrected 
for the detector efficiency and acceptance loss. In 
this section, we discuss separately the single-electron 
efficiencies and electron pair efficiencies. 

1. Single-electron efficiency 

The single-electron efficiency is determined by the 
product of the TPC tracking efficiency £xpc, the 
TOF matching efficiency £xof, and the electron 
identification efficiency £ e ro: 

£ e = £TPC X £xOF X £ e ID (8) 

The TPC tracking efficiency, £xpc, was evaluated 
via the standard STAR embedding technique. In the 
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unLikeS(mix)/LikeS(mix) LikeS(same)/unLikeS(mix) unLikeS(mix)/LikeS(mix) LikeS(same)/unLikeS(mix) 



FIG. 19. Ratios of unlike-sign to like-sign mixed-event distributions (denoted as the acceptance-difference correction 
factor) and acceptance corrected like-sign background to mixed-event unlike-sign background distributions for different 
centralities and pt regions. 



butions of signal pairs compared to the inclusive unlike- 
sign (open symbols) and reconstructed background pairs 
(black line) in 200 GeV Au+Au minimum-bias colli¬ 
sions. The insert shows the signals of the tf> and u> vector 
mesons. Panel (b): Signal-to-background ratios in p+p 
and Au+Au collisions. 


embedding process, simulated electron tracks with a 
certain phase space definition were generated and 
then passed through the STAR detector geometry 
for the 2010 (2011) configuration using the Geant 
model. Next, the simulated detector signals were 
mixed with real data to have a realistic detector oc¬ 
cupancy environment. The mixed signals were pro¬ 
cessed with the same offline reconstruction software 
that was used for the real data production. The 
tracking efficiency was studied by comparing the re¬ 
constructed tracks with the simulated input tracks. 
The input number of simulated tracks (5% of to¬ 
tal event multiplicity) were constrained to prevent a 
sizable impact on the final single-track efficiency. 

The electron track TOF-match efficiency, £tof, 
was obtained from real data samples. Due to the 
limited pure electron statistics, we first used a pure 
pion sample in order to deduce the TOF-match ef¬ 
ficiency. Pure pion samples were selected based on 
a TPC dE/dx cut. We assume the TOF-match ef¬ 
ficiencies for different particle species are similar in 
the pt region where dE/dx cannot distinguish dif¬ 
ferent particle species. Pure electron samples were 
selected to cross-check the efficiency scale differences 
between electrons and pions due to the decay loss of 
pions between the TPC and the TOF detectors as 
well as other effects. Electrons (or positrons) from 
photon conversion or 7r° Dalitz decays were iden¬ 
tified by invariant mass and topological techniques 
and used as the high purity samples. 

The TPC tracking and TOF matching efficien¬ 
cies were calculated differentially in three dimensions 
(pr,r),<J)). The pion TOF matching efficiency was 
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also calculated in (px > P, < /0 while a same scaling fac¬ 
tor, which accounts for the TOF matching efficiency 
difference between pions and electrons, was used for 
all ( 77 , <f>) bins due to limited statistics. The choice 
of the binning in ( 77 , <f>) dimensions shows a negligi¬ 
ble effect in the -integrated final dielectron pair 
efficiency. 

The electron identification cut efficiency, £ e iD, in¬ 
cludes two components: efficiency due to the TOF 
1//3 cut (ep) and efficiencies due to the dE/dx PID 
selection criteria (EdEdxPro)- 

E e ID = ep X EdEdxPID , . 

(9) 

EdEdxPID — E n dEdx X E na - e 

Pure electron samples were used to study the TOF 
1//3 distributions. In order to estimate the 1//3 effi¬ 
ciency, Ep, we applied two methods to the 1//3 distri¬ 
butions: a realistic function fit and direct counting. 
The difference in the results from the two methods 
was included in the systematic uncertainty. 

The dE/dx PID selection efficiency, EdEdxPID, in¬ 
cludes the efficiency due to the cut on both the num¬ 
ber of dE/dx points and na e which is used to select 
the electron candidates. The cut efficiency on the 
number of dE/dx points, e n dEdx, was deduced using 
the pure pion samples in the real data. The results 
from the electron sample were consistent with those 
from pions in the region allowed by the statistics of 
the samples used. Then the efficiency from the pion 
samples was used in the final efficiency calculation in 
three dimensions (px,p, <(>)■ The na e cut efficiency, 
£ncr e) was deduced via the same steps as described 
in Section III-D for calculating the electron purity 
and hadron contamination. With the extracted na e 
gaussian mean position and width values, the PID 
cut efficiency was calculated under the selection cri¬ 
teria described in Section III-C. 

In the upper left panel of Fig. [2U Expc(pt), 
£tof (pt), £ndEd x (pt) and their product are shown 
for tracks in minimum-bias collisions. These ef¬ 
ficiencies are averaged over \rj\ < 1 and 2-7T in az¬ 
imuth. The ratios of £tpc x £tof x E n dEdx at differ¬ 
ent centralities are shown in the bottom left panel of 
Fig- HU The ep, £ nrTe and their product as a function 
of momentum are shown in the upper right panel of 
Fig. HU The centrality dependence of ep x e n(Te is 
shown in the bottom right panel. 


2. Electron pair efficiency 

The dielectron pair efficiency was evaluated from 
the single-electron efficiency in the following two 
ways: 


• Toy Monte Carlo simulation, which used the 
virtual photons as the input and let them de¬ 
cay into dielectrons isotropically. 

• Cocktail simulation, which used the hadronic 
cocktail (see Secttion G) as input including the 
correlated heavy-flavor decay electrons from 
Pythia simulations [ 32 J. 

In the final dielectron spectra, we have experimen¬ 
tal ambiguities in separating heavy-flavor decayed 
dielectron yields from medium-produced dielectron 
yields (including contributions from both hadronic 
and partonic sources). Furthermore, the heavy- 
flavor decay dielectron production is not known in 
heavy-ion collisions due to possible medium modifi¬ 
cations of the heavy-flavor correlations when com¬ 
pared to those in p+p collisions. We used these 
two methods to estimate our dielectron pair effi¬ 
ciency. The single-electron efficiencies, described in 
the previous section, were folded in for each daughter 
track in a full three dimensional (px,p, <(>) momen¬ 
tum space. The pair efficiency and acceptance was 
finally calculated in (M ee ,px). 

Shown in Fig.[22]are the dielectron pair efficiencies 
in the STAR acceptance (pf- >0.2 GeV/c, \rj e \ < 1) 
with \i/ee\ < 1- The difference in pair efficiency in 
the STAR acceptance between these two methods is 
small, ranging from about 3% at low px, down to 
about 1% at high px- And due to statistical lim¬ 
its of the cocktail simulation for the dielectron from 
heavy-flavor decay, we use the pair efficiency calcu¬ 
lated from the virtual photon decay in this analysis 
and include the difference between these two meth¬ 
ods in the systematic uncertainty. 

The (jy\ pair cut efficiency was evaluated using a 
7T° embedding sample in which simulated 7r° parti¬ 
cles with enriched Dalitz decays were embedded into 
the real data. The efficiency was calculated after 
re-weighting the input 7r° yield with a realistic px 
distribution (details in the next part). We also used 
a pure virtual photon decay convoluted with the de¬ 
tector resolution for this calculation. The difference 
was included as the systematic uncertainty of the </>v 
pair cut efficiency. 

In Fig. [231 the e + e _ pair efficiencies are shown 
as a function of pair px in different mass regions. 
In the high px/mass region the efficiency is almost 
constant as the single track efficiency turns stable at 
high px (see Fig. HU)- The px-integrated e + e pair 
efficiencies as a function of pair mass within STAR 
acceptance in Au+Au collisions at ^/snn = 200 GeV 
are shown in Fig. 1241 The pair efficiency without the 
4>v cut is also plotted, which contributes only in the 
very low mass region. 
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Transverse momentum (GeV/c) momentum (GeV/c) 


FIG. 21. The single-electron efficiency as a function of pr or p in the pseudo-rapidity range of |p| < 1 in Au+Au 
collisions at ^/snn = 200 GeV. In the left panel solid points represent the combined TPC tracking, TOF matching 
and TPC ndEdx efficiencies. In the right panel the solid points represent the combined TPC na e cut and TOF 
1//3 cut efficiencies. Gray bands indicate the systematic uncertainties. The bottom panels show these efficiencies in 
various centrality bins compared to the minimum-bias data. 



FIG. 22. e + e~ pair efficiency as a function of pair 
mass in different px regions calculated from two different 
methods. 
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G. Hadronic cocktails 

Dielectrons as measured by the detector originate 
from all stages in the evolution of heavy-ion colli¬ 
sions. These pairs include the decay products of 
long-lived particles which typically decay after they 
have frozen out of the medium. The contributions 


FIG. 23. e + e _ pair efficiency as a function of pair pr for 
different mass regions. The dashed lines represent ±10% 
difference from the unit, the solid lines show a constant 
fit to the data. 


























































20 



FIG. 24. (Color online) The px-integrated e + e“ pair ef¬ 
ficiency as a function of pair mass within the STAR ac¬ 
ceptance in Au+Au minimum-bias (black) and central 
(red) collisions at ^/snn = 200 GeV. The open markers 
represent the corresponding efficiencies before the pho¬ 
ton conversion rejection. The gray bands depict the sys¬ 
tematic uncertainties. 


in the final dielectron spectrum can be evaluated as 
long as their yields at freeze-out are known. 

The simulation process for constructing the con¬ 
tributions from such decays in Au+Au collisions, of¬ 
ten referred to as the hadronic cocktail, is similar to 
what has been done in p+p collisions and reported 
in j3lj. The cocktail simulations only contain the 
hadron form-factor decays in the vacuum at freeze- 
out. Cocktails included in our calculation contain 
contributions from decays of 7 r°, rj , ?/, w, (j>, J/ip, 4 '> 
cc, bb as well as from Drell-Yan (DY) production. A 
vacuum p meson calculation is included separately 
when discussing the data compared to cocktail with 
the vacuum p. For the hadron decay calculations, 
the input rapidity distributions are assumed to be 
flat within |y| < 1. The input yields dN/dy within 
this rapidity window as well as the px distributions 
are discussed below. 

The charged pion yields at 200 GeV Au+Au colli¬ 
sions have been accurately measured in the STAR 
acceptance [H, Is3 |. The input 7 r° spectrum is 
taken as the averaged yield between STAR’s 7 r + and 
n~ measurements. Other light hadron yields in¬ 
clude the p meson, measured by PHENIX for py > 
2 GeV/c |35], and meson data from STAR (3(|. 
These hadron spectra together with hadron spectra 
(K^, K[ j, and A) measured by STAR and PHENIX 
were simultaneously fit to a core-corona based Tsal- 
lis Blast-Wave (TBW) model [37:j where the core de¬ 
scribes the Au+Au bulk production and the corona 
describes the hard scattering contribution from p+p 
like collisions. 


In Fig. [55j the simultaneous fit results for all in¬ 
put hadron spectra are shown except for J/ip- The 
J/4 contribution is not considered as a component 
of the bulk medium. The cocktail input for J/ip was 
taken from the measurement by the PHENIX collab¬ 
oration [38} . For light hadrons, the TBW functions 
provide good parameterizations to these measured 
spectra. For those hadron cocktail components with¬ 
out corresponding direct measurements (e.g. low px 
p, 77 ', oj), we use the same core TBW parameters ob¬ 
tained from the fit and predict the spectral shapes 
for each of these unknown components, shown as 
solid curves in Fig. 1251 The low px 77 spectrum was 
fixed by requiring the match with the measured data 
points at px> 2GeV/c, while the dN/dy of 77 ' me¬ 
son was taken with the same values as used in the 
PHENIX publication [Tt}. The same set of TBW 
parameters from the simultaneous fit were used to 
generate the ui spectrum and the dN/dy was tuned 
to match our dielectron yield in the w peak region. 

Additional corrections were applied to account for 
the differences in centrality and rapidity windows 
between the input hadron spectra and our dielec¬ 
tron measurements. The measured pion yields were 
calculated in the rapidity window of \y\ < 0.1 in 
Ref. [Hj and \y\ < 0.5 in Ref. [H}. We used the 
pion rapidity distributions from the Hijing calcu¬ 
lations and scaled the measured pion yields down 
by 3% to obtain the px spectrum in the rapidity 
window of \y\ < 1. This correction factor was also 
included in the uncertainty of the input 7 r° dN/dy. 
The different centrality windows matter when tak¬ 
ing the minimum-bias data from PHENIX measure¬ 
ments, done in 0-92% centrality, and compare those 
to our results which are for 0-80% centrality. We 
corrected for this difference using the measured 7 r° 
dN/dy values in 0-92% and 0-80% centralities by the 
PHENIX experiment [39] . 

The correlated charm, bottom and Drell-Yan 
contributions were obtained from Pythia calcula¬ 
tions [32] and scaled by the number of binary col¬ 
lisions in Au+Au collisions for the default cocktail 
calculations. We used Pythia version 6.419 with 
parameter settings: MSEL=1, PARP(91) (( kj _)) = 
l.OGeV/c and PARP(67) (parton shower level) = 
1.0. This setting was tuned to match our measured 
charmed meson spectrum in p+p collisions [40| . 
The input charm-pair production cross section per 
nucleon-nucleon collision was also taken from charm 
meson measurements (13, EH- We used the same 
PYTHIA setting to calculate the dielectron yields 
from correlated bottom decays and from the Drell- 
Yan production. The input bottom and Drell-Yan 
production cross sections are: a/// = 3.7 pb, = 
42 nb. 

The p meson contribution is expected to be mod- 
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FIG. 25. (Color online) Invariant yields of mesons in 
Au+Au collisions at ^/snn = 200 GeV. The solid lines 
show the simultaneous TBW fit to the measured data 
points and the TBW predictions for r/, rf , and t o with 
the same set of fit parameters. The dashed lines show 
the same parametrization to the measured J/ip spectrum 
and the predicted ip' spectrum as in ll7 ]. 


ified due to a strong coupling to the hot QCD 
medium created in heavy-ion collisions. Therefore, 
the p meson was not included in our default cocktail 
calculations. In the comparison between our mea¬ 
sured dielectron spectra and the cocktail calculations 
including the vacuum p , we used the p meson mea¬ 
surements in peripheral collisions by STAR [42j and 
assumed a similar p/ir ratio in order to extrapolate 
to other centrality selections. The mass spectrum 
of the vacuum p —> e + e~ was taken the same line 
shape as reported in our dielectron measurement in 
p+p collisions [31 ] . 

Table QH] summarizes all sources of the hadron 
cocktail and their decay branching ratios. The 
TBW d?} parametrizations were used to describe 
the input hadron px distributions, shown in Fig. 1251 
The resulting e + e _ pair mass distributions from the 
individual sources are normalized by the respective 
decay branching ratios and measured yields dN/dy. 
Additional scaling parameters for various centrality 
bins are listed in Appendix [A] 

The mass spectra reported in this paper are not 
corrected for the STAR detector resolution. It is 
very challenging to precisely reproduce the momen¬ 
tum resolution in the STAR TPC simulation pack¬ 
age in the high luminosity RHIC environment due to 
various distortion effects in the TPC detector. In¬ 
stead a data-driven method was used to obtain the 
diclectron mass line shape in the cocktail simulation. 

Based on the full detector simulation, the recon¬ 
structed electron p™ c probability distribution at a 
given input pJ(? c was parametrized with a double 


Crystal Ball function [46|, defined as: 


P(Pt C ,Pt C ) °c 


Ax {B — R)~ n , R<-a 

e - , — a < R < /3(10) 

Cx (D + R )~ m , R> /? 


with 
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( 11 ) 


where n = 1.29, a = 1.75, m = 2.92, and pS = 
1.84. The value of p = —0.001 is slightly shifted 
from 0 due to the electron energy loss in the detector 
material as STAR tracking accounts for the energy 
loss assuming all tracks are pions. 

The pt resolution is taken as <r PT /pt and assumed 
to follow the form: 



(axp T ) 2 + 
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/3 


P_ 
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Pt 

yjptp + TO 2 
( 12 ) 


For electrons /3 ~ 1. 

We used the J/ip signal which has the most statis¬ 
tics and tuned the parameters a and b in the Eq. [T5] 
to get the best match to the J/ip signal distribu¬ 
tion. The two parameters were found to be a = 
6.0 x 10“ 3 c/GeV and b = 8.3 x 10" 3 . 


H. Systematic uncertainties 

The major sources of systematic uncertainty that 
contribute to the final result in this analysis include: 

1. Normalization factor for mixed-event distribu¬ 
tions 

2. Residual correlated background 

3. Like-sign/unlike-sign acceptance difference 
correction 

4. Hadron contamination 

5. Efficiency and acceptance corrections 
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TABLE III. Input yields of various cocktail components for 0-80% minimum bias Au+Au collisions at ^/snn 
200 GeV. 


source 

B.R. 

dN/dy or a 

Uncertainty 

Reference 

7r u —>■ qee 

1.174 x 10~ 2 

98.5 

8% 

STAR [33, 341 

r] —> 7 ee 

7 xlO” 3 

7.86 

30% 

PHENIX [17, 351 

v' -»• 

'yee 

4.7 xlO -4 

2.31 

100% 

PHENIX [171, STAR [311 

H 

ee 

4.72 xlO -5 

16.7 

42% 

STAR [42J 

UJ —>■ 

ee 

7.28 xlO -5 




oj — > 7r°ee 

7.7 xlO -4 

9.87 

33% 

STAR [43J 

cp ->• 

ee 

2.95 xlO -4 




r]ee 

1.15 xlO -4 

2.43 

10% 

STAR [36] 

JH>- 

-»■ ee 

5.94 xlO -2 

2.33 xlO -3 

15% 

PHENIX [38] 

ipi — 

ee 

7.72 xlO -3 

3.38 xl0“ 4 

27% 

PHENIX [44, 45] 

CC —» 

ee 

1.03 xlO” 1 

da cc /dy = 171pb 

15% 

STAR [41] 

bb —» ee 
DY ee 

1.08 xlO -1 
3.36 xlO -2 

Cpp = 3.7 pb 
&PP = 42 nb 

30% 

30% 

Pythia[32] 

Pythia[32] 


The systematic uncertainty of the background of 
diclectron pairs was further separated in two mass 
regions, where we chose different background sub¬ 
traction methods (see Eq. 0. 

In the mass region of M > 1.0GeV/c 2 , we ob¬ 
tained the signal by subtracting the mixed-event 
unlike-sign background plus the residual correlated 
background. The normalization of the combinato¬ 
rial background, applied to the mixed-event unlike- 
sign distribution, is determined by comparing the 
like-sign same-event and mixed-event distributions. 
The statistics of the total like-sign pair in the nor¬ 
malization region is the dominant systematic uncer¬ 
tainty. We also chose different normalization ranges 
varying between the mass range of 1.2— 2.0 GeV/c 2 . 
Other sources that we considered include the nor¬ 
malization method and the slight asymmetry be¬ 
tween the total number of mixed-event unlike-sign 
and like-sign pairs. For the normalization method, 
we chose a different method compared to what was 
described in Section III-E, in this way we normal¬ 
ize the mixed-event unlike-sign distribution to the 
acceptance-corrected same-event like-sign distribu¬ 
tion. Table ED summarizes the contributions for 
each of the individual components to the system¬ 
atic uncertainty of the normalization factors in min¬ 
imum bias as well as for various other centralities 
from 200 GeV Au+Au collisions. 

The uncertainty in the residual correlated back¬ 
ground was already mentioned in Section III-E.5. In 
the data-driven approach, the statistical uncertainty 
in determining the ratio of like-sign and mixed-event 
unlike-sign r(M ee ,pT) was used as the systematic 
uncertainty. The contribution to the final dielectron 
mass spectrum in minimum-bias collisions is about 
10% from 1 GeV/c 2 to 3 GeV/c 2 . 

In the low mass region, M ee <1.0 GeV/c 2 , we ob¬ 
tained the signal by subtracting the acceptance cor¬ 


rected like-sign background, in which the acceptance 
difference correction between like-sign and unlike- 
sign pairs was calculated using mixed-event distri¬ 
butions. Different event mixing methods by vary¬ 
ing the different event categories and event pool 
sizes were chosen, and the largest deviations between 
these methods are used in the uncertainty calcu¬ 
lation. The acceptance correction done in the 2D 
( M ee , px) plane may suffer from limited statistics. 
The difference between the results calculated using 
the 2D acceptance correction and using the ID ( M ee 
only) acceptance correction was included in the sys¬ 
tematic uncertainty as well. 

The electron candidates contain a small amount 
of hadron contamination, which may be correlated 
with other particles ( e.g. from resonance decays) 
and thus contribute to the final signal spectrum. 
To estimate this contribution, we first selected pure 
pion, kaon and proton samples with stringent TOF 
to 2 limits. We randomly picked hadrons from these 
pure samples according to the estimated hadron con¬ 
tamination levels in both the total amount and the 
Pt differential yields, creating a hadron contamina¬ 
tion candidate pool. The analysis procedure used in 
the dielectron analysis was applied to that pool to 
estimated the e — h and h — h correlated contribu¬ 
tions. 

The estimated hadron comtamination evaluated 
from e — h and h — h correlated contributions com¬ 
pared to the dielectron signal is shown in Fig. 01] 
Overall, the relative contribution to the final spec¬ 
trum is <5% between 1 GeV/c 2 and 3GeV/c 2 . 

The systematic uncertainties on the raw dielec¬ 
tron invariant-mass spectra for minimum-bias colli¬ 
sions are summarized in Fig. 1271 As a conservative 
estimation, we took the sum of each individual com¬ 
ponent as the total systematic uncertainty. 
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TABLE IV. Systematic uncertainties on normalization factors of mixed-event distributions for minimum-bias collisions 
and various centralities. The total number of e + e“ pairs in minimum-bias collisions is ~ 3.7 x 10 7 , and for central 
collisions is ~ 7.0 x 10 7 (2010 data). 



Like-sign pairs 

Choice of N.R. 

Norm, method 

LS/US pair difference 

Total 

MinBias 

4.9 x 10 -4 

2.1 x 10 -4 

1.0 x 10 -4 

2.4 x lO -3 

0.05% 

0-10 % 

3.4 x 10" 4 

1.4 x 10~ 4 

5.6 x 10 -5 

1.7 x 10~ 5 

0.04% 

10-40% 

6.6 x 10“ 4 

3.2 x 10" 4 

1.2 x 10“ 4 

3.1 x 10“ 5 

0.07% 

40-80% 

2.2 x 10" 3 

5.2 x 10~ 3 

5.2 x 10" 4 

9.8 x 10 -5 

0.56% 



FIG. 26. Estimated hadron contamination from e — h 
and h—h contributions due to the finite contaminated 
hadrons in the elecron sample compare to the e + e“ pair 
signal in 200 GeV Au+Au collisions. 



FIG. 27. (Color online) Systematic uncertainties 
of raw dielectron invariant-mass spectrum in Au+Au 
minimum-bias collisions from various contributing 
sources. The direct sum of each individual component 
was taken as the total systematic uncertainty, shown as 
the solid curve. 


For the reported dielectron yields in the STAR 
acceptance, the systematic uncertainty due to the 
efficiency correction includes uncertainties on the 
single-track efficiency, the pair efficiency evaluation 
method, and the pair cut (cf>v) efficiency. Table fVl 
summarizes each individual component and their 
contributions to the total uncertainty of the single- 
track efficiency. The individual component contri¬ 
butions were determined by varying track selection 
cuts and comparing distributions between data and 
MC for the uncertainty on the TPC tracking effi¬ 
ciency (nHitsFits, DCA, etc.). The uncertainties on 
the TOF matching efficiency, the TOF PID cut ef¬ 
ficiency, and the ndEdxFits cut efficiency were eval¬ 
uated by comparing the results obtained from the 
pure electron samples from photon conversion and 
7T° Dalitz decay. The difference between a realistic 
function fit and direct counting methods of the TOF 
1//3 distribution was also included in the uncertainty 
of the TOF PID cut efficiency. The electron pair effi¬ 
ciency evaluated from single tracks was described in 
Section III-F.2. Due to the unknown relative contri¬ 
butions between the correlated charm decays and the 
medium contribution, two extreme calculations were 
used as conservative estimates for the systematic un¬ 
certainty. This uncertainty is mostly constrained to 
the intermediate and high mass regions of the mass 
spectrum and ranges from about 3% at low px down 
to 1% at high px- The systematic uncertainty of the 
(j)y pair cut efficiency was evaluated by taking the 
difference between the calculations from the 7r° em¬ 
bedding sample and the virtual photon decay sam¬ 
ple, which is about 3% at M ee < 0.05GeV/c 2 . The 
systematic uncertainty of the pair efficiency due to 
different methods is 5%. Finally the total systematic 
uncertainty of the electron pair efficiency is ~ 13%. 


IV. RESULTS AND DISCUSSION 

A. Dielectron mass spectrum in minimum-bias 
collisions 

The dielectron yields measured in the STAR ac¬ 
ceptance (p|. > 0.2GeV/c, |p e | < 1, and \y ee \ < 1) 
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TABLE V. Systematic uncertainties on single-track effi¬ 
ciency. 



component 

uncertainty 


nHitsFits 

4.0% 

TPC 

DCA 

ndEdxFits 

2.5% 

1.0% 


n CF e 

2.0% 

TOF 

matching 

i/d 

1.0% 

3.0% 

Total 


6.1% 


have been corrected for the dielectron reconstruction 
efficiencies. The efficiency correction was done in 
Px an d M ee . The px-integrated efficiency-corrected 
dielectron mass spectrum dN/dM at midrapidity 
| p ee | < 1 in the STAR acceptance from 0-80% 
Au+Au minimum collisions at ^/snn = 200 GeV is 
shown in Fig. (2U The data are compared to the 
hadronic cocktail simulations without (upper left 
panel) and with (upper middle and upper right pan¬ 
els) the vacuum p contribution. The vertical bars 
on the data points depict the statistical uncertainty, 
while the green boxes represent the systematic un¬ 
certainty. The ratios of the data over the cocktail 
simulations are shown in each of the bottom pan¬ 
els. The yellow band around unity indicates the 
uncertainties on the cocktail calculations. Those 
are mainly determined by the uncertainties on the 
dN/dy and the decay branching ratios for each of 
the individual sources. 

A few more remarks about cocktail calculations 
are in order: 

• Since the p mesons are strongly coupled to the 
medium in Au+Au collisions, their contribu¬ 
tion is considered part of the medium dilepton 
emission and depends on the properties of the 
medium. We only included the vacuum p con¬ 
tribution as a reference here. In the default 
hadronic cocktail calculations, the p contribu¬ 
tion is omitted in order to allow for possible 
in-medium p contributions depicted by model 
calculations. 

• Correlated charm contributions included in 
the cocktail are the number-of-binary colli¬ 
sions (IVbin) scaled p+p results calculated from 
Pythia. 

• Other hadron contributions are described in 
Section III-G. 

Comparing the measured data points to the 
hadronic cocktail calculations in the LMR, an 
enhancement can be observed in the mass region 
between 0.30 and 0.76GeV/c 2 . This enhance¬ 
ment cannot be fully explained by the expected 


vacuum p meson contribution as shown in the 
right plot of Fig. [28] The data, integrated in 
the mass region of 0.30—0.76 GeV/c 2 , is a factor 
of 1.76 ± 0.06 (stat.) ± 0.26 (sys.) ± 0.29 (cocktail) 
larger than the model cocktail without the vacuum 
p contribution. This enhancement factor is signifi¬ 
cantly lower than what has been reported from the 
dielectron measurement in the PHENIX detector 
acceptance 0- 

Detailed comparisons of the differences between 
the STAR and PHENIX experimental acceptances 
and cocktail simulations are unable to account for 
the measured enhancement difference. These details 
are described in Appendices |B] and [Cj 

In the intermediate mass region (IMR), the cock¬ 
tail simulations are dominated by correlated charm 
pair decays which are calculated from Pythia simu¬ 
lations. The simulations generally describe the data 
but run slightly below the data points, allowing for 
additional source contributions. The uncertainty on 
the charm production cross section da cc /dy at mid¬ 
rapidity, which is used for the normalization of this 
contribution, is around 15%. More precise measure¬ 
ments in this mass region of both the total charm 
cross-section as well as the correlation in Au+Au col¬ 
lisions are needed to either verify or rule out signifi¬ 
cant contributions from other sources, such as QGP 
thermal radiation. 


B. Comparison to models 

One major motivation for measuring dileptons 
is the study of chiral symmetry properties of the 
QCD medium created in the heavy-ion collisions. 
Restoration of the spontaneously broken chiral sym¬ 
metry will lead to modification of the vector meson 
(short-lived p meson in particular) spectral func¬ 
tions, which are accessible via dilepton measure¬ 
ments. There are two chiral symmetry restora¬ 
tion scenarios commonly used in calculations: (a)the 
drop of the pole mass or degeneracy of vector and 
axial-vector mesons due to the reduced ( qq) conden¬ 
sate 0; and (b) broadening of the spectral func¬ 
tion due to many body collisions in the vector-meson 
dominance [§, [ill, S3]. Both scenarios will intro¬ 
duce an enhancement in the mass region below the p 
mass compared to the spectral function in vacuum. 
Precision measurement from the NA60 experiment 
demonstrated that the broadened p scenario can re¬ 
produce the low mass dilepton enhancement data at 
SPS energy [ju|, while the dropping mass scenario 
failed to describe the data. It is anticipated that the 
hadronic medium at top RHIC energy is similar to 
that created at SPS energy, thus the dilepton pro¬ 
duction in the LMR region are comparable between 
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FIG. 28. (Color online) Invariant mass spectrum in the STAR acceptance (pf. > 0.2GeV/c, |p e | < 1, and \y ee \ < 1) 
from y /snn = 200 GeV Au+Au minimum-bias collisions. The mass spectrum is compared to the hadronic cocktail 
simulations without (upper left panel) and with (upper middle) the vacuum p contribution (upper right is an expanded 
version of upper middle below M ee of 1.4 GeV/c 2 ). The vertical bars on data points depict the statistical uncertainties, 
while the systematic uncertainties are shown as green boxes. Yellow bands in the bottom panels depict the systematic 
uncertainties on the cocktail. The dashed line indicates the charm decay dielectron contribution from Pythia [32] 
calculations and scaled with Abm- 


SPS and RHIC. 

The QGP contribution to the dilepton spectra 
has often been calculated perturbatively via Born 
q + q annihilation at leading order. Various ap¬ 
proaches have been studied to take into account 
high-order contributions at finite T — ps HU]. The 
QGP contribution is expected to become sizable for 
M > 1.5 GeV/c 2 at top RHIC energies due to a well 
established partonic phase. 

There have been many model calculations for di¬ 
electron production at RHIC, with particular focus 
on the low mass region. We group these models into 
two categories and describe their features and pre¬ 
dictions separately below. 

Category I: Macroscopic effective many-body the¬ 
ory models. In these models, the dilepton produc¬ 
tion in the hadronic medium is calculated via elec¬ 
tromagnetic correlators based on the Vector-Meson 
Dominance Model (VDM) approach. Assuming 
a thermal equilibrated hadronic medium, dilepton 
rates are determined by the p meson propagator in 
the medium, which depends on the interactions of 
the p with mesons and baryons in this medium at 
finite T and ps- It has been shown that the re¬ 
sulting broadened p spectral function is mostly due 
to the interactions with the baryons rather than 
the mesons [49l45l| . Thus, the medium total-baryon 
density, and not the net-baryon density or pg, is the 
critical factor in determining the dielectron yield in 
the heavy-ion collisions at these energies. 

Dilepton production in the partonic phase is 
mostly calculated via perturbative qq annihilation 


with some improved corrections. It has been demon¬ 
strated in these calculations that the dilepton rates 
from the hadronic medium, extrapolated bottom-up 
to T c , should be equivalent to the rates from the par¬ 
tonic medium, extrapolated top-down to T c . This 
is referred to as the “parton-hadron” duality HU- 
The final resulting diclectron yields for observation 
are calculated via the integral over the full space- 
time evolution for this medium. We have chosen 
one model calculation from Rapp [52| from this cate¬ 
gory in the following comparisons to our data. Some 
of the key ingredients in this model calculation are 
listed below: 

• Initial spectral functions were fixed using the 
measurements from e + e _ collision data 11]. 

• Space evolution was chosen to be a cylindrical 
expanding fireball [HI] 

• The latest Lattice QCD Equation-of-State 
(EOS) was used, in particular lower T c (170 
MeV) and T c h (160 MeV) values were cho¬ 
sen in the calculations shown here, which are 
slightly different compared to previous calcu¬ 
lations from this same model. 

• QGP radiation from the partonic phase was 
updated as well, and using the choice of the 
latest Lattice QCD EOS. 

There are several other model calculations avail¬ 
able in this category: some models chose differ¬ 
ent initial spectral functions HU, and several of 
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them used the space-time evolution obtained from 
either ideal or viscous hydrodynamic model calcula¬ 
tions [50|, .51] . Calculations from these models show 
similar results compared to Rapp’s model and pro¬ 
vide reasonable descriptions of the low-mass excess 
observed in our diclectron data in 200 GeV Au+Au 
minimum-bias collisions. 

Category II: Microscopic transport dynamic mod¬ 
els. We chose the Parton-Hadron String Dynamic 
(PHSD) transport model from this category when 
comparing to our data in the following sections. 
The PHSD transport approach incorporates the rel¬ 
evant off-mass-shcll dynamics of the vector mesons 
and an explicit partonic phase in the early hot and 
dense reaction region as well as the dynamics of 
hadronization frit . It allows for a microscopic study 
of the various dilepton production channels in non¬ 
equilibrium matter. In the hadronic sector, PHSD is 
equivalent to the HSD transport approach that has 
been used for the description of p +A and A+A colli¬ 
sions from SIS to RHIC energies. It reproduces fairly 
well the measured hadron yields, rapidity distribu¬ 
tions, and transverse momentum spectra [54| . The 
dilepton radiation by the constituents of the strongly 
interacting QGP is produced via: (i) basic Born 
q + q annihilation, (ii) gluon Compton scattering 
(q/q + g —> 7* + q/q), and (iii) quark/anti-quark an¬ 
nihilation with the gluon bremsstrahlung in the final 
state (q + q —> g + 7*) . Dilepton production in these 
partonic channels is calculated with off-mass-shell 
partons using a phenomenological parametrization 
for the quark and gluon propagators in the QGP. 

The PHSD model has been used to calculate 
the dielectron yields in the STAR acceptance and 
it shows a fair agreement with our preliminary 
data [55l] . 

Detailed comparisons of the model calculations 
with the data are shown in Fig. [55] In the LMR, 
the data and model calculations are in a fairly good 
agreement. In the IMR, the charm contribution is 
the most important component of e + e _ spectrum. 
We discuss the effect of possible modification of this 
component in Section IV-F. 


C. pt dependence 

To gain more insight into dielectron production, 
we studied the px dependence of the dielectron yields 
in comparison to the hadron cocktail and model cal¬ 
culations. In different px regions, comparisons to 
hadron cocktails require precise knowledge of the 
light hadron production in a wide px region. De¬ 
tails of the cocktail calculations on the px shape of 
input particle are described in Section III-G. 

The measured diclectron yields within STAR ac¬ 


ceptance in each individual px region as well as 
the total expected hadron cocktail contributions are 
shown in the left panel of Fig. [3D] Note that the 
correlated charm contributions, which become very 
important in the mass region from 0.5—3.0GeV/c 2 , 
were all taken from the IVbin scaled Pythia calcu¬ 
lations. The ratios of data over cocktail calculations 
as a function of M ee for several transverse momen¬ 
tum ranges are shown in the right panels of Fig. 1301 
For comparison, the theoretical model calculations 
in each px window are included as well. The en¬ 
hancement factor with respect to the hadronic cock¬ 
tail does not change significantly in these px bins. 
Both theoretical models are able to reasonably de¬ 
scribe the LMR excess in all px bins. 

We quantify the px dependence by comparing the 
measured diclectron yields with the cocktail in each 
mass window within the STAR acceptance, the re¬ 
sults from Au+Au 0-80% minimum-bias collisions 
at 200 GeV are shown in Fig. [3T] The left panel 
shows the measured data points (markers) together 
with cocktail calculations (dashed lines). The ratios 
of the data over the cocktail are shown in the right 
panels for different mass windows. The data points 
and the cocktail calculations are in good agreement 
throughout the measured px range up to 2 GeV/c 
in the mass regions of the 7 r° (up to 0.15 GeV/c 2 ), 
the ui/(j) (0.76—1.05GeV/c 2 ), and the J/V' mesons 
(2.8—3.5 GeV/c 2 ). In the LMR region, particularly 
in the mass region of 0.30—0.76GeV/c 2 , we see that 
the relative enhancement in the data compared to 
the cocktail has no significant px dependence. Ta¬ 
ble [Vi] summarizes the enhancement factors for each 
px bin. In the IMR region, cocktail calculation can 
describe the data reasonably well. Due to the large 
uncertainty on the correlated charm contribution 
there is little constraint on other possible dilepton 
contributions, e.g., QGP thermal radiation. 

TABLE VI. ThepT dependence of dielectron yields, mea¬ 
sured in the STAR acceptance, and the enhancement fac¬ 
tor with respect to the hadronic cocktail in the invariant 
mass region of 0.30—0.76 GeV/c 2 . 


p T (GeV/c) 

yield (xlO d ) 

yield/cocktail 

0 -0.5 

115 + 0.09 + 0.20 

1.71+0.12 + 0.29 

0.5-1.0 

158 + 0.07 + 0.27 

1.56 + 0.07 + 0.27 

1.0-1.5 

066 + 0.03 + 0.11 

1.81+0.09 + 0.29 

1.5-2.0 

024 ± 0.02 ± 0.04 

2.65 + 0.16 + 0.44 


D. Centrality dependence 

The dielectron spectra are studied in various cen¬ 
trality bins (0-10%, 10-40% and 40-80% ). The 
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FIG. 29. (Color online) Dielectron mass spectrum in 200 GeV minimum-bias Au+Au collisions compared to the 
hadron cocktail plus the hadronic medium and partonic QGP contributions calculated from Rapp (upper left panel) 
and PHSD (upper right panel) models. Yellow bands in the bottom panels depict systematic uncertainties on the 
cocktail. 
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FIG. 30. (Color online) Left panel: Invariant mass spectra from ^/snn = 200 GeV Au+Au minimum-bias collisions 
in different pr ranges. The solid curves represent the cocktail of hadronic sources and include the charm-decayed 
dielectron contribution calculated by Pythia scaled by A+; n . Right panels: The ratio of dielectron yield over cocktail 
for different pr bins, and the comparison with model calculations. The gray boxes show the systematic uncertainties 
of the data. Yellow bands depict systematic uncertainties on the cocktail. 
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FIG. 31. (Color online) Left panel: The integrated dielectron yield as a function of pair pr in different invariant mass 
ranges compared with cocktail. The solid lines represent the yield of the cocktail in different mass ranges , while the 
gray bands show the systematic uncertainties of the data. Right panels: The ratio of dielectron yield over cocktail for 
different mass ranges as a function of pair pr . The yellow bands show the systematic uncertainties of the cocktail. 
The gray bands show the systematic uncertainties of the data. 


left panel of Fig. [35] shows the dielectron spectra 
in these centrality bins compared to cocktail calcu¬ 
lations. The ratios of the data to the cocktail are 
presented in the right panels. Model calculations are 
also included in the right plots. In Fig. [33] we quan¬ 
tify the measured yields as a function of centrality 
by means of 7V part for different mass windows. 

In the LMR, particularly in the mass region 
0.30—0.76 GeV/c 2 , the observed enhancement factor 
of the dielectron yield with respect to the cocktail 
does not show a significant centrality dependence 
within current uncertainty. Both theoretical models 
can reasonably reproduce the centrality dependence 
of this observed enhancement in the LMR. Table fVTTI 
summarizes the enhancement factors for each cen¬ 
trality bin. 

In Fig. [M] we overlay the dielectron mass spectra 
from minimum-bias and the most central (0-10%) 
collisions for which we are able to achieve sufficient 
statistics for direct comparisons. The 7V par t-scaled 
spectra are plotted in the upper panel, and the ratio 
between them is plotted in the bottom panel. The 
measured ratio is consistent with unity in the 7r° in¬ 
variant mass region, indicating that the production 
scales with N palt . The ratio starts to increase in 
the mass around 0.5—1.0GeV/c 2 . This observation 
is consistent with a picture in which the correlated 


charm contribution starts to be a dominant source 
in this mass region while charm quark production at 
RHIC is expected to rather scale with N\,i n . Addi¬ 
tionally, in this invariant mass range the in-medium 
p meson contribution from the hadronic medium is 
expected to increase faster than 7V part when moving 
towards central collisions based on model calcula¬ 
tions [52]. In the IMR, the data indicate there is 
potentially a systematic change in the mass spectra 
when comparing the minimum-bias and central colli¬ 
sions. This is suggestive of a possible modification of 
charmed hadron production or other contributions 
such as thermal radiation. To quantify the differ¬ 
ence, exponential fits were performed to the mass 
spectra in central and minimum-bias collisions and 
the resulting exponential slopes differ by ~1.5cr. 

TABLE VII. The centrality dependence of dielectron 
yields, measured in the STAR acceptance, and the en¬ 
hancement factor with respect to the hadronic cocktail 
in the invariant mass region of 0.30—0.76 GeV/c 2 . 


centrality 

yield (x 10 3 ) 

yield/cocktail 

0 - 10% 

13.63 ± 1.01 ± 2.06 

2.03 ±0.15 ±0.31 

10 - 40% 

4.81 ±0.22 ±0.71 

1.63 ±0.08 ±0.24 

40 - 80% 

0.85 ±0.03 ±0.12 

1.51 ±0.06 ±0.22 

0 - 80% 

3.87 ±0.13 ±0.57 

1.76 ±0.06 ±0.26 
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FIG. 32. (Color online) Left panel: Invariant mass spectra from ^/snn = 200 GeV Au+Au collisions in different cen¬ 
tralities. The solid curves represent the cocktail of hadronic sources and include charm decay dielectron contribution 
from Pythia scaled by A^m. Right panel: The ratio of dielectron yield over cocktail for different centrality. The gray 
boxes show the systematic uncertainties of the data. Yellow bands depict systematic uncertainties on the cocktail. 
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FIG. 33. (Color online) Left panel: The integrated dielectron yield in different centralities compared with the cocktail. 
The solid lines represent the yield of cocktail in different invariant mass ranges. The gray bands show the systematic 
uncertainties of the data. Right panel: The ratio of dielectron yield over cocktail for different centralities, and the 
comparison with model calculations. The yellow bands show the systematic uncertainties of the cocktail. The gray 
bands show the systematic uncertainties of the data. 
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FIG. 35. (Color Online) Mass spectrum of the ex¬ 
cess (data - cocktail) in the low-mass region in Au+Au 
minimum-bias collisions compared to model calculations. 
Green brackets depict the total systematic uncertainties 
including those from cocktails. 


FIG. 34. (Color Online) Upper panel: Dielectron in¬ 
variant mass spectra from minimum-bias and the most 
central (0-10%) collisions that we are able to achieve 
most statistics at present. For direct comparison, the 
spectra are scaled with the number of participant nucle¬ 
ons (Ap ar t). The solid line represents the cocktail for 
minimum-bias collisions. Bottom panel: The ratio of 
the Ap ar t-scaled dielectron yield between minimum-bias 
and the most central collisions. The gray boxes show the 
systematic uncertainties of the data. 


E. Low mass excess yields 

We subtracted the cocktail contribution from the 
measured diclecton mass spectrum to obtain the di¬ 
rect excess yields, shown in Fig. [35] for the mass re¬ 
gion of 0.3-1.4 GeV/c 2 . The cocktail simulations 
used in the subtraction include the correlated charm 
contributions from Pythia assuming the Abi n scal¬ 
ing. A possible charm de-correlation leads to a negli¬ 
gible modification of the cocktail spectra in the mass 
region around 0.5 GeV/c 2 as shown in Fig. [34] The 
obtained excess spectra in Au+Au minimum-bias 
collisions are compared to model calculations in Fig. 

M 

The systematic uncertainty across all the data 
points are highly correlated. We utilized the modi¬ 
fied y 2 -test [56J to quantify the comparison between 
the data and the model calculations; the results are 
summarized in Table I yTh! The vacuum p plus QGP 
scenario in Rapp’s implementation cannot describe 
our data well. The calculations, including the broad¬ 
ened p-meson scenario plus QGP contribution from 
both Rapp and PHSD, have reasonable agreements 
with our data. 

Next, we studied the centrality dependence of the 


excess yields. In Fig.[36l the integrated excess yields 
scaled by A part as a function of centrality (A part ) are 
shown in the p-like mass region (0.30-0.76 GeV/c 2 ). 
In the same figure, the w-like (0.76-0.80 GeV/c 2 ) 
and 0-likc (0.98-1.05 GeV/c 2 ) dielectron yields are 
plotted. For both sets, the yields were scaled by 
A par t and the cocktail subtraction was not applied 
in this range. The w-like and the 0-like dielectron 
yields show an A part scaling while the p-like dielec¬ 
tron excess yields increase faster than A part as a 
function of centrality. The dashed curve depicts a 
power fit (oc A part ) to the p-like dielectron yields 
with the cocktail subtracted. The fit result shows 
a = 0.44 ± 0.10 (stat.+uncorrelated sys.), indicat¬ 
ing the dielectron yields in the p-like region are sen¬ 
sitive to the QCD medium dynamics, as expected 
from p medium modifications in theoretical calcula¬ 
tions m, 57]. 

TABLE VIII. Reduced y 2 for model calculations com¬ 
pared to the excess data in the invariant mass region of 
0.3-1.0 GeV/c 2 . 


Model 

X7ndf 

p- value 

Rapp 

vacuum p + QGP 

41.3/8 

2.4x 10 Y 

Rapp 

broadened p + QGP 

8.0/8 

0.32 

PHSD 

broadened p + QGP 

16.5/8 

0.040 


F. Correlated charm contributions 

The correlated charm contributions start to play 
an important role in our measured dielectron yields 
above 0.5 GeV/c 2 and dominate the cocktail in the 
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FIG. 36. (Color Online) The yields scaled by A par t for 
the p-like region with the cocktail subtracted, and the cj- 
like and the 0-like regions without cocktail subtraction 
as a function of -/V par t. Model calculations are included as 
red solid and dot-dashed curves, while the dashed curve 
depicts a power-law fit to the yield/lV par t for the p-like 
region with the cocktail subtracted. Systematic uncer¬ 
tainties from the data are shown as grey boxes, and the 
green brackets depict the total systematic uncertainties 
including those from cocktails. For clarity, the cu-like and 
0-like data points are slightly horizontally displaced. 


intermediate mass region. So far, no measurement of 
charm correlation in the low transverse moment re¬ 
gion at RHIC exists in either p+p or A+A collisions. 
Single-charm hadron spectra or their decay (“non¬ 
photonic’’) electron spectra have been measured in 
p+p 0, [EH an d Au+Au collisions [U, @ • We re¬ 
lied on the Pythia model to create the correlated 
charm pairs and then calculate the decay-electron 
pair distributions. 

In p+p collisions, with a tuned Pythia set¬ 
ting: MSEL=1, PARP(91) {(k±)) = 1.0 GeV/c and 
PARP(67) (parton shower level) = 1.0, we have 
shown that this can reproduce the measured sin¬ 
gle Z+meson px spectrum from 0.6 — 6GeV/c jioj . 
The dielectron mass spectrum calculated with this 
Pythia tune also showed a good agreement with 
our measurement in the IMR in p+p collisions at 
200 GeV. However, the limited statistics in p+p col¬ 
lisions do not allow us to determine whether Pythia 
can produce the correct D — D correlation. 

In Au+Au collisions, we have observed that high- 
px electrons are strongly suppressed compared to 
p+p collisions. In the low px region, various model 
calculations indicate that the single-charm spectrum 
can be modified due to interactions between charm 
quarks and the hot and dense medium [60]. Con¬ 
sequently, the D—D correlation seen in p+p colli¬ 
sions will be modified, or even be completely washed 
out H2- To study their impact on the dielectron 
spectrum, we chose the following different config¬ 


urations for the charm px spectra and correlation 
functions to construct D—D pairs. 

(a) Keep the direct Pythia calculation which was 
used in our default cocktail calculations. 

(b) Keep the momentum magnitude of charm de¬ 
cay electrons in Pythia, but randomly select 
the azimuthal angle direction. In this case, 
the angular correlation between two electrons 
is completely washed out. 

(c) Randomly sample two electrons with the single 
electron px, 77, 4> distributions from Pythia 
calculation. In this case, the correlation be¬ 
tween the two electrons is completely washed 
out. 

(d) Based on (c), but sample the px of each elec¬ 
tron track according to the modified px dis¬ 
tribution based on the non-photonic electron 
Raa measurement in Au+Au collisions (59j . 
The electron Raa(pt) was parametrized us¬ 
ing the following function, with px in units of 
GeV/c. 

Raa(j>t) = 4 63 + e ( PT -o. 62 )/i.o 6 ~ °- 22 ' ( 13 ) 

All these calculations were scaled with iVbin in 
each centrality bin to obtain the correlated charm 
mass spectra. The correlated charm mass spectra 
for the above four cases in the most central (0-10%) 
Au+Au collisions, and a comparison with the mea¬ 
sured data, are shown in Fig. [37] The total cocktail 
shown is still calculated based on the default Pythia 
correlations. The figure shows that both the modi¬ 
fication in electron momentum and the smearing in 
azimuthal angular correlation make the dielectron 
mass distribution steeper. Calculations for case (d) 
seem to be closer to the data points in the mass re¬ 
gion of 1—3 GeV/c 2 , thus indicating a possibly mod¬ 
ification of charmed hadron production in central 
Au+Au collisions that is worthy of further experi¬ 
mental investigations. We also calculated the slope 
parameter in the transverse mass spectrum for each 
of the aforementioned cases, as is shown in Fig. l38l 

G. Low mass vector meson yields 

The low mass vector meson (ui and (j>) yields have 
been extracted from the diclectron decay channel 
through this analysis. The results reported here are 
from combined data taken in RHIC year 2010 and 
2011 runs. The measured <f> yields are consistent 
with the results from a recent STAR publication @1 • 
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FIG. 37. (Color Online) Correlated charm contributions 
to the dielectron mass spectra for different assumptions 
of the correlation strength. The simulations are com¬ 
pared to the measured dielectron spectrum in the most 
central (0-10%) Au+Au collisions. The total cocktail 
shown is calculated using the default Pythia correla¬ 
tions. The insert plot shows the comparsion between 
the cocktail sums with above four different charm con- 
tribtuions and the measured spectrum in the IMR. 



FIG. 38. (Color Online) Slope parameters T e g versus in¬ 
variant mass for dielectrons from charm hadron decays. 
Different lines show the outcome from Pythia calcula¬ 
tions assuming different levels of correlations between 
charm pairs. 


Figure [321 shows the invariant mass distributions of 
the vector mesons uj and <f> from ^/snn = 200 GeV 
Au+Au minimum-bias collisions. The signal spec¬ 
tra are reconstructed by subtracting the normalized 
mixed-event unlike-sign background (Section III-E3) 
from the inclusive same-event unlike-sign e + e _ dis¬ 
tribution. A Breit-Wigner function plus a second 
order polynomial function are used to fit the invari¬ 
ant mass distributions. The second order polynomial 
function is used to describe the residual background. 
In addition, we use the vector meson uj{(j)) invariant 
mass distributions (line-shapes) directly from cock¬ 
tail simulations (Breit-Wigner plus Gaussian func¬ 
tions) to fit the signal. As described in Section III- 
G, the detector momentum resolution in the cock¬ 
tail simulation was estimated by tuning the simula¬ 
tion to match the J/ip signals in the data. The line 
shapes from this tuned simulation for the uj and <f> 
mesons reproduce the signal well. The difference be¬ 
tween these two methods is included in the system¬ 
atic uncertainty of the raw yield. Figure [40] shows 
the uj and </> invariant mass distributions in different 
Pt regions. 

Although the mass and width of vector mesons 
could be modified due to interactions with the hot 
and dense medium, the observed uj and </> spectra 
from the detector will have little sensitivity to such 
an effect. The lifetimes of uj and <f> mesons are much 
longer than the typical lifetime of the medium cre¬ 
ated in high energy heavy-ion collisions. Therefore, 
the freeze-out uj and (f> mesons will dominate the 
observed yields. We obtained the widths and mass 
positions of uj and (f> signals from data and compared 
them to the values from the PDG as well as from our 
simulations, shown in Fig. 1411 

The mass positions of the uj and cf> mesons from 
our data generally agree with the PDG values, with 
a slight shift towards lower values. This is mainly 
because the STAR tracking algorithms account for 
the energy loss assuming pions. The observed mass 
shift (1-2 MeV/c 2 ) for uj and 4> mesons are within 
the uncertainties of the particle energy loss correc¬ 
tion in our GEANT simulations. The widths of the 
mass distribution are larger than the PDG values 
as expected, due to detector resolution effects. A 
tuned simulation, using the J/i/j mass distribution 
(described in Section III-G), can reproduce well the 
observed signal widths for uj and (f> mesons in the 
full px region reported here. Because of uncertainty 
in the description of materials, including accessory 
components in the STAR detector system, we in¬ 
cluded the difference between the tuned simulation 
and the GEANT simulation in the width calculation 
as part of our systematic uncertainty. Since the mass 
and width are well reproduced by the tuned simu¬ 
lation, we fixed the mass and width with the value 









33 



FIG. 39. (Color Online) ui and <f> meson invariant mass distribution from ^/snn = 200 GeV Au+Au minimum-bias 
collisions after subtraction of the combinatorial background using the mixed-event method. The blue and red lines 
depict two functions used for the signals in the fit. A second-order polynomial function is used to describe the residual 
correlated background. 
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FIG. 40. (Color Online) pt dependence of the uj and meson invariant mass distributions from ^/snn = 200 GeV 
Au+Au minimum-bias collisions. 


from the tuned simulations when using the Breit- 
Wigner function to extract the pp differential yield. 

In order to present the final p^-differential in¬ 
variant cross section, the raw vector meson yields 
are corrected for the detector’s acceptance and ef¬ 
ficiency. Figure 02] shows the total detector accep¬ 
tance and efficiency as a function of pp for ui —> 
e + + e~ and (f> -+ e + + e~. In Fig. 03] the final 
Pt differential invariant yields are shown for ui and 
</> from = 200 GeV Au+Au collisions at mid¬ 

rapidity (| y | < 1). The systematic uncertainties in¬ 
clude the detector efficiency uncertainty and the raw 
signal extraction uncertainty. The latter is derived 
from changing the fit range, the function used for 
describing the background and the method used to 
extract the yields. The <f> spectrum measured from 


e + e _ decays is consistent with the previous results 
measured from the hadronic decay channel {4> —> 
K + + K ~) [36]. Also included in the figure are the 
TBW-model filTt fit to the previous (f> -+ K + + K~ 
data points as well as a prediction of the ui spectrum 
with the same set of parameters obtained from the 
simultaneous fit to all available light hadrons (see 
Section III-G). The TBW prediction describes the 
measured w spectrum well. The measured dN/dy 
for the w meson is 8.46±0.67(stat.)±1.59(syst.), and 
for the (f> meson is 2.20±0.10(stat.)±0.34(syst.). 
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FIG. 41. (Color Online) The widths and mass positions 
of the cj and <j> signals from data compared to the values 
from the PDG and the full Geant simulation. Boxes 
on the data points depict the systematic uncertainties. 
Gray bands represent the uncertainty of the simulation. 



FIG. 42. (Color Online) The efficiency and acceptance 
correction factor as function of px for mid-rapidity ui and 
<t> mesons. 



FIG. 43. (Color Online) The pr distributions of the to 
and 4> meson invariant yields from ^/snn = 200 GeV 
Au+Au minimum-bias collisions. The solid black and 
blue lines are from the Tsallis Blast-Wave functions as 
described in Section fill Gl 


V. SUMMARY 

We have reported STAR measurements of di¬ 
electron yields at mid-rapidity in Au+Au collisions 
at -y/sNN = 200 GeV. The measured dielectron 
yields within the STAR acceptance (defined by 
Pt > 0.2 GeV/c, \p e \ < 1, and \y ee \ < 1) show an 
enhancement when compared to hadronic cocktail 
calculations in the mass region below M^. The 
enhancement factor, integrated over the mass region 
of 0.30 — 0.76 GeV/c 2 and the full px acceptance, 
is 1.76 ± 0.06 (stat.) ± 0.26 (sys.) ± 0.29 (cocktail). 
Further systematic measurements show that this 
enhancement factor has a mild centrality and px 
dependence. A vacuum p spectral function cannot 
fully describe the measured dielectron mass spec¬ 
trum in this mass region. This enhancement factor 
is significantly lower than what has been reported 
by PHENIX. We have compared the STAR and 
PHENIX cocktail simulations and applied PHENIX 
azimuthal acceptance. We found that neither differ¬ 
ences in the acceptance nor the cocktail simulations 
can explain the difference in the enhancement factor 
measured by the two experiments. 

We compared our results to model calculations in¬ 
cluding an effective many-body model (Rapp) and a 
microscopic transport model (PHSD). Both models 
invoked an in-medium modified p spectral function 
through the interactions with mesons and baryons 
in the bulk medium. Both can reproduce the low- 
mass excess in our data reasonably well, including 
the px and centrality dependences. A power-law fit 
to the excess yield in the p-like region as a function of 
Apart gives a power of 1.44+0.10. We noted that the 
many-body model calculations have successfully ex- 
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plained the SPS low-mass dilepton data. These find¬ 
ings could indicate that the property of the hadronic 
medium that governs the dilcpton production in the 
low-mass region is similar at top SPS and top RHIC 
energies despite the difference in center-of-mass en¬ 
ergies of more than an order of magnitude. Di¬ 
electron measurements from the RHIC beam energy 
scan program will offer a unique opportunity to fill 
the energy gap between the SPS and RHIC and sys¬ 
tematically evaluate the energy dependence of di¬ 
electron production. 

We also reported the measurement of w meson 
production, and <j> meson production through the di¬ 
electron decay channel in Au+Au collisions at ^/snn 
= 200 GeV. The observed signal widths and mass 
positions are well reproduced in Monte Carlo sim¬ 
ulations. The measured <j) invariant yield spectrum 
through the e + e - decay channel is consistent with 
the previously published STAR measurement based 
on the K + K~ decay channel. The to invariant yield 
spectrum can be well reproduced by Tsallis Blast- 
Wave model predictions which use the same set of 
parameters obtained from a simultaneous fit to all 
other available light hadrons. This indicates that 
the ui mesons behave much like the bulk medium, 
with similar radial flow velocity. 

The understanding of the dielectron production in 
the mass region of 1—3GeV/c 2 is currently limited 
both statistically and systematically. We reported 
the inclusive dielectron yields which include the con¬ 
tribution from correlated charm decays. However, 
at this time we do not know the characteristics of 
the charm contribution in Au+Au collisions. The 
diclectron data from 0-80% minimum-bias collisions 
can be fairly well described by the number-of-binary- 
collisions scaled p+p contribution based on Pythia 
calculations. The ratio between the central and 
minimum-bias spectra in the mass region of 1-3 
GeV/c 2 shows a moderate deviation from the Win 
scaling (1.8er deviation for the data point at 1.8-2.8 
GeV/c 2 ). This could be indicative of the modifica¬ 
tion of the correlated charm contribution or the exis¬ 
tence of other contributing sources in Au+Au colli¬ 
sions. The difference in the mass region 1-3 GeV/c 2 , 
if confirmed in future measurements with a better 
precision, would constrain the magnitude of the de- 
correlating effect on charm pairs while traversing 
the QCD medium and/or other possible dielectron 
sources, e.g., QGP thermal radiation, from central 
Au+Au collisions at RHIC. 
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Appendix A: Centrality-dependent Cocktail 
Simulation Inputs 


When comparing to the measured spectra, the 
hadron cocktails were simulated for each of the 
corresponding centrality bins (0-10%, 10-40% and 
40-80%). The centrality dependence of the input 
hadron px distributions were obtained from the sim¬ 
ilar Tsallis Blast-Wave function fit to the available 
data, including <j>, etc. For other hadrons with 
no available measurements, we use the Tsallis Blast- 
Wave predictions for the input px distributions. The 
input dN/dy for all the components in each cen¬ 
trality bin were then scaled with the relative pion 
yields, R n , with respect to minimum-bias collisions. 
Correlated-charm contributions are scaled by the 
number of binary collisions for a given centrality. All 
of these scale factors are summarized in Table IIXI 


TABLE IX. Scale factors for centrality dependent cock¬ 
tails 


centrality 

dN(-7r)/dy 

Rtt 

(AW) 

0-80% 

98.49 

1 

291.90 ± 20.46 

0-10% 

279.2 

2.834 

941.24 ± 26.27 

10-40% 

131.1 

1.331 

391.36 ±30.21 

40-80% 

30.45 

0.309 

56.62 ± 13.62 
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FIG. 44. (Color online) Single electron/positron track 
density distributions using the STAR, data selected 
within the PHENIX azimuthal angle acceptance. 


Appendix B: Detector Acceptance Effect 


The STAR mid-rapidity detectors cover the full 
azimuth (0< cj> < 2n) for \r]\ <1 while the PHENIX 
central arms (used for the dielectron analysis) cover 
about 2 x 7 t/ 2 at \rj\ < 0.35. To investigate the im¬ 
pact of the detector acceptance effect on the final di¬ 
electron mass spectrum, we tried to narrow down the 
single-track acceptance cut to match the PHENIX 
acceptance as best as possible. We acknowledge that 
fully reproducing another experiment’s acceptance 
is virtually impossible due to subtle differences in 
detector structures and performances. Instead, the 
STAR data is selected with the PHENIX azimuthal 
angle acceptance cut. Due to the limited statistics, 
we cannot further reduce the pseudorapidity window 
to match the respective PHENIX range. In addition, 
we also expect the physics is not significantly differ¬ 
ent between \rj\ < 0.35 and |fy| < 1 rapidity ranges 
for 200 GeV collisions. 

As a result of the magnetic field, the signal track 
</> acceptance varies with px ■ We used the kinematic 
acceptance cut presented in the PHENIX publica¬ 
tion H3 : 


‘/’min < 4> + 9 


‘/’min < </> + q 


fcpc 

Pt 

fcRICH 

Pt 


— i/max 
A cj )max 


(Bl) 


where kuc = 0.206 rad-GeV/c and fcRicH = 
0.309 rad-GeV/c represent the effective azimuthal 
bend to DC and RICH, respectively. One arm cov¬ 
ers the region from </> mi n = -^7r to </> max = ^7r 
, while the other arm covers from </> m ; n = ^tt to 

, _ 19 

rmax — Jg • 

The electron candidate occupancy distributions 
for STAR data selected with the PHENIX (f> accep¬ 
tance cut are shown in Fig. 1441 The upper panel 
shows the regular <j) versus px for negative charged 
particles, while the bottom panel shows 1/px ver¬ 
sus <j) for both charges. The plots show that while 
we can capture the basic acceptance structure, the 
inner fine structure within this azimuthal angle ac¬ 
ceptance may be slightly different due to different 
detector structure for both experiments. 

With the electron candidates selected, we then 
carried out the same analysis procedure as described 
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FIG. 45. Ratios of pair distributions for electron can¬ 
didates selected with the PHENIX cf> acceptance. Pan¬ 
els (a), (b), (c): ratios of like-sign distributions between 
same event and mixed-event. 


in Section III. In Fig. IT51 panels (a), (b), (c), the ra¬ 
tio is shown for like-sign distributions between same 
events and mixed events from which we determine 
the normalization factor of the mixed-event unlike- 
sign distribution for the combinatorial background. 

We compared the acceptance difference correc¬ 
tion factor between the results with and without the 
PHENIX (j) acceptance, as shown in Fig. [46] One 
can clearly see that the cj> acceptance cut changes 
the pair acceptance between like-sign and unlike- 
sign pairs significantly in the low-mass region, and 
the maximum of this ratio appears around M ee ~ 
0.5 GeV/c 2 . 

The combinatorial background is subtracted from 


the inclusive unlike-sign pair distribution to obtain 
the raw signal, then the raw signal distribution is 
corrected for the detector efficiency. Finally we ob¬ 
tained the signal dielectron invariant mass spectrum 
from 200 GeV minimum-bias Au+Au collisions and 
compared it to hadronic cocktail simulations, shown 
in the left panel in Fig. [17] In the low-mass region of 
0.30-0.76 GeV/c 2 , we observed an enhancement of a 
factor of 2.4±0.37(stat.)±0.38(sys.)±0.29(cocktail) 
when comparing the measured yield to the hadronic 
cocktail. Selecting our data within the PHENIX <f> 
acceptance does not appear to reproduce the large 
enhancement factor in the low-mass region observed 
by the PHENIX collaboration [13] ■ 

We added in the medium dielectron contributions 
from theoretical model calculations. The right plot 
of Fig. [TO shows the data compared to the cock¬ 
tail plus the broadened p spectrum in the hadronic 
medium and QGP thermal radiation. The partic¬ 
ular calculation that is included here is only valid 
for M < 1.5GeV/c 2 . The medium contribution 
from this model (hadronic p and QGP radiation) 
describes the observed low-mass excess very well. 


Appendix C: Cocktail comparison between 
STAR and PHENIX 

In this appendix, we compare the cocktail sim¬ 
ulation results between PHENIX and STAR. The 
details of the light hadron decays and Dalitz de¬ 
cays into dielectrons are described in [31]]. We used 
the p+p input yields for cocktail calculations in this 
comparison and we folded in the PHENIX accep¬ 
tance filter, described in Eq. m Next, we com¬ 
pared the output to the PHENIX p+p cocktail cal¬ 
culations. 

The comparison for the total cocktail summed 
yield is shown in Fig. 1481 The comparisons for each 
cocktail component are shown in Fig. [T)] We see 
that the cocktail calculations from both experiments 
agree reasonably well. There are some small differ¬ 
ences in the rp to, (f> Dalitz decay distributions which 
can be attributed to different choices of decay form 
factors in these Dalitz decays. 
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FIG. 46. (Color online) Left panel: unlike-sign/like-sign pair acceptance difference correction factor with the PHENIX 
4> acceptance (black stars) compared with the full acceptance (red circles). Right panel: dielectron pair mass distri¬ 
butions of the STAR data within the PHENIX <f> acceptance; the inclusive unlike-sign distribution (black), the same 
event like-sign distribution (blue), and the mixed-event unlike-sign (red), like-sign (magenta) distributions. 
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FIG. 47. (Color online) Left panels: efficiency corrected invariant mass spectra (blue solid dots) calculated using 
the STAR data filtered with the PHENIX azimuthal angle acceptance. The data points are compared to cocktail 
simulations shown as curves and the lower left panel shows the ratio of data to the cocktail sum. Right panels: the 
same data points compared to theoretical model calculations of medium vector meson and QGP contributions from 
Ref. [52j. The bottom right panel shows the ratio of data to the sum of the cocktail and the theory calculations of 
medium contributions. 
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(black dashed) for p+p 200 GeV collisions. 
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